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aiACK  PROPAGATION  AND  FABRIC  CONTROL 


ON  THE  STATIC  AND  DYNAMIC  BEHAVIOR  OF  COHESIVE  SOILS 


by 

Adel  S.  Saada  and  Gary  F.  Bianchini 


RESEARCH  OBJECTIVES  AND  ACCOMPLISHMENTS 

To  conduct  parallel  studies  on  strength  behavior,  kinematics  of  deformation  and  soil 
fabric,  in  saturated  clay  soils;  with  particular  emphasis  being  placed  on  crack  and  damage 
propagation,  bifurcation,  and  localization  in  the  form  of  shear  bands.  The  following 
questions  were  asked  in  the  original  research  proposal: 

1.  How  does  the  fabric  and  the  strength  of  a  cross  anisotropic  clay  change  with  the 
inclination  of  the  principal  stresses  on  the  axis  of  symmetry  of  the  material? 

2.  What  is  the  effect  of  a  crack  or  notch  and  of  anisotropy  on  the  dynamic  modulii  E 
and  G  and  on  the  damping  ratio  ,  at  small  strains  in  fast  cyclic  tests? 

3.  How  does  the  fabric  change  in  the  vicinity  of  a  crack  propagating  in  Mode  n  during 
fast  and  slow  cyclic  testing? 

4.  How  does  a  damage  zone  propagate  under  static  loading  as  compared  to  its 
propagation  under  fast  loading  and  cyclic  loading?  What  is  the  influence  of  the 
amplitude  and  of  the  number  of  cycles? 

5.  What  is  the  nature  of  the  densification  strips  in  the  damage  zone?  What  happens  to 
the  clay  clusters  in  those  strips? 

6.  How  do  cracks  and  damage  zones  coalesce  in  case  of  more  than  one  crack 
propagating  initially  in  Mode  II? 

7.  How  is  the  fabric  modified  as  the  damage  zones  propagate  and  coalesce? 

All  the  questions  were  answered.  In  addition,  the  observations  and  measurements 
made  along  the  shear  bands  during  their  propagation  provide  the  necessary  ingredients  for 
the  development  of  the  appropriate  constitutive  equations.  Testing  with  hollow  cylindrical 
specimens  allows  one  to  apply  a  whole  spectrum  of  stress  states  to  a  material.  TTiis  is  the 
first  time  that  initiation  and  propagation  of  shear  bands  were  measured  in  a  material 
subjected  to  a  tri-dimensional  system  of  stress.  The  fact  that  we  are  dealing  with  a 
relatively  weak  material  allows  us  to  obtain  excellent  results  while  using  simple  equipment. 


STUDENTS  PARTICIPATION 


Four  students  were  involved  of  whom  three  left  after  spending  relatively  short 
periods  of  time  working  on  this  research.  Mr.  Liqun  Liang  the  fourth  student  joined  the 
project  a  year  and  a  half  ago  and  contributed  substantially  to  the  conduct  of  the 
experiments.  The  bulk  of  the  work,  however  was  done  by  the  principal  investigators.  Mr. 
Liang  will  be  using  the  results  to  develop  a  finite  elements  code,  most  probably  based  on 
elasto-plasticity,  to  model  the  bifurcation  and  shear  banding  phenomena  in  his  Ph.D.  thesis. 

PUBLICATIONS  &  PRESENTATIONS 

Two  papers  are  being  prepared  for  publication.  This  is  a  project  in  which  one  has  to 
wait  to  the  end  in  order  to  reach  definite  conclusions. 

A  presentation  will  be  made  in  Paris  in  June  of  1992  during  a  Franco-American 
Workshop  on  Strain  Localization  and  Interfaces;  the  workshop  being  supported  and 
sponsored  by  the  National  Science  Foundation  and  the  French  "Centre  National  pour  la 
Recherche  Scientifique"  (CNRS). 

In  addition,  a  presentation  will  be  made  at  the  November  1992  ASCE  New  York 
Convention  during  a  session  on  Fracture  Mechanics  Applied  to  Soils. 
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ABSTRACT 


Discontinuities  in  the  form  cracks  or  fissures  and  inclusions  are  often  present  in 
natural  clays.  They  serve  as  stress  concentrators  when  loads  are  applied  to  the  material. 
Such  concentrations  result  in  the  cracks  advancing,  often  surrounded  and  preceded  by  a 
propagating  damage  zone.  As  the  crack  propagates,  the  damage  may  be  in  the  form  of  one 
or  more  shear  bands,  which  may  play  the  part  of  new  stress  concentrators  and  blimt  the 
action  of  the  original  crack.  This  report  examines  some  of  the  phenomena  associated  with 
the  presence  of  cracks  in  overconsolidated  clays  as  well  as  in  the  less  brittle  normally 
consolidated  clays.  Differences  between  isotropic  and  anisotropic  materials,  and  the  level 
at  which  serious  modifications  take  place  in  the  fabric  of  the  material  are  noted. 
Observations  made  with  a  surface  analyzer  as  well  as  with  a  transmission  electron 
microscope  indicate  that  the  changes  in  fabric  mostly  occur  at  the  level  of  the  cluster  or 
flock.  The  influence  of  the  cracks  and  of  the  shear  bands  on  the  kinematics  and  strength  of 
the  test  specimens  is  studied  and  discussed. 
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CHAPTER  1 


INTRODUCTION  AND  SHORT  REVIEW  OF  PAST  WORK 


The  mechanical  response  of  clays  is  dictated  by  their  fabric  as  well  as  by  flaws  or 
discontinuities  that  may  exist  in  the  material.~These  flaws,  often  in  the  form  of  cracks  (or 
notches),  result  in  shear  bands  that  spread  from  the  crack  and  form  a  damaged  zone  in 
which  most  of  the  deformation  is  localized.  Bifurcation  theories  have  been  and  are  being 
used  to  predict  the  occurrence  of  shear  bands  based  on  the  pre-localization  constitutive 
relations  of  the  material. 

Many  branches  of  mechanics  must  be  brought  together  to  analyze  the  phenomena  at 
hand.  The  following  briefly  summarizes,  or  better  lists,  key  investigations  conducted  by 
various  researchers;  the  hope  being  that  at  some  time  a  coherent  formulation  will  be  found 
to  express  the  various  observed  phenomena. 

1.1  In  the  Area  of  Shear  Band  Observation: 

Morgenstem  and  Tchalenko  (1967)  used  thin  sections  of  carbowax  impregnated  clay 
to  study  the  behavior  of  a  clay  tested  in  a  direct  shear  device.  They  defined  two  types  of 
discontinuities,  namely  strain  and  displacement,  and  bring  out  the  kinematic  restraint 
imposed  by  the  configuration.  They  also  adopt  a  sequential  numbering  system  in  referring 
to  the  shear  bands.  Their  detailed  examination  of  the  development  of  kink  bands  and 
compression  textures  lead  to  the  notion  of  basal  plane  slip  as  the  basic  mechanism  in  shear 
induced  structures. 
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Houlsby  and  Wroth  (1980)  classify  kineniatic  discontinuities  and  show  that  a 
velocity  discontinuity  does  not  necessarily  have  to  be  a  line  of  zero  extension;  whether  it  is 
or  not  depends  on  the  adopted  constitutive  relation. 

Balasubramaniam  (1976)  working  with  clays  used  X-rays  to  study  the  pattern  of 
deformation.  He  found  the  strains  to  be  uniform  up  to  75  percent  of  the  peak  value. 

Arthur  and  Dunstan  (1982)  used  radiography  to  show  that  in  granular  non-cohesive 
material  rupture  layers  formed  gradually  by  local  layers  of  higher  dilation  in  pre-failure 
strain  align  along  pre-failure  incremental  non-extension  directions,  axes  of  stress  and  strain 
rate  being  generally  coincident.  Rupture  layers,  formed  at  failure  in  homogenous  samples 
of  soil  sheared  slowly  under  homogeneous  stresses  with  fixed  principal  stress  directions,  are 
oriented  at  a  specific  angle  according  to  soil  particle  size  between  the  limits  of  (x/4  - 
072)  and  (jr/4  -  P/2);  the  limits  being  determined  by  the  plane  strain  values  of  0  ’  and 
the  angle  of  dilation. 

Scarpelli  and  Wood  (1982)  used  radiography  to  study  the  orientation  of  the  shear 
bands  in  direct  shear  tests  conducted  on  sand.  They  suggest  that  the  degree  of  constraint 
perceived  by  the  sand  will  influence  the  particular  bifurcation  that  it  adopts  at  any 
particular  location.  Where  the  sand  sees  freedom  it  may  adopt  the  Vardoulakis-Arthur 
solutions.  Where  the  imposed  constraint  is  larger,  as  in  the  direct  shear  device  the  shear 
band  would  aim  for  the  Roscoe  solution  and  follow  directions  of  zero  extension. 

Desrues,  J.  (1983)  used  stereophotogrametry  to  study  localization  in  granular 
materials.  For  sands  Desrues  accurately  measured  the  direction  of  shear  bands  due  to 
plane  strain  compression  of  sand  samples.  He  provided  local  measurements  of  the 
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dilatancy  angle  and  showed  that,  under  certain  circumstances,  the  local  dilatancy  angle 
within  the  shear  band  was  twice  as  large  as  the  global  dilatancy  angle  calculated  from  the 
total  volume  change  measurements. 

12  In  the  Area  of  Bifurcation  Theories ; 

Vardoulakis  (1982)  in  his  review  of  the  theoretical  and  experimental  work  done  in 
the  area  of  stability  and  bifurcation  of  soil  samples,  states  that  research  in  this  area  has 
been  primarily  influenced  by  the  work  of  Hill  (see  Hill  and  Hutchinson,  1975).  Hill  and 
Hutchinson  (1975)  based  their  work  on  Biot’s  pioneering  studies  (1965)  and  Ariaratnam 
and  Dube/s  (1969)  studies  of  difruse  bifurcation  modes  under  plane  strain  conditions  for 
elasto-plastic  incompressible  materials.  They  studied  in  great  detail  bifurcation  of 
incompressible  materials  in  the  plane-tension  tests  covering  a  broad  spectrum  of 
constitutive  assumptions.  They  covered  both  diffuse  and  localized  shearing  modes.  Plane 
compression  tests  were  analyzed  by  Young  (1976),  and  more  recently  Needleman  (1979) 
extended  this  work  to  pressure  sensitive  incompressible  materials. 

Rudnicki  and  Rice  (1975)  studied  localization  by  comparing  an  isotropic  hardening 
model  with  a  yield  vertex  model  for  a  pressure  sensitive,  dilatant  material.  They  showed 
that  the  yield  vertex  model  provides  a  better  description  and  that  non  normality  permits 
non  uniqueness  with  positive  hardening.  Both  these  features  are  shown  to  be  destabilizing 
and  to  strongly  influence  the  resulting  predictions  for  localization  by  comparison  to 
predictions  based  on  classical  plasticity  idealizations,  involving  normality  and  smooth  yield 
surfaces. 
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Rice  (1976)  discusses  the  localization  of  plastic  deformation  into  a  shear  band  as  an 
instability  of  plastic  flow  and  a  precursor  to  rupture.  The  material  is  assumed  rate- 
independent  and  critical  conditions  are  sought  for  its  constitutive  relations  to  allow  a 
bifurcation  from  a  smoothly  varying  deformation  into  a  highly-concentrated  shear  band 
mode.  He  essentially  explores  a  bifurcation  approach  in  which  concentrated  deformation 
occurs  simultaneously  (or  near  simultaneously)  at  all  points  of  the  zone  of  localization. 
The  results  are  applied  to  rigid-plastic,  elastic-plastic  and  plastically  dilatant  with  pressure 
sensitive  yielding  materials. 

Rice  and  Rudnicki  (1980)  make  a  distinction  between  continuous  and  discontinuous 
bifurcation.  In  the  first  further  plastic  deformation  is  assumed  to  occur  Uo  h  inside  and 
outside  the  zone  of  localization;  while  in  the  second  elastic  unloading  occurs  outside  the 
zone  of  localization  while  continued  elastic-plastic  deformation  occurs  within  it. 
Discontinuous  bifurcations  become  possible  at  the  same  time  continuous  ones  occur.  The 
elastic  unloading  will  cause  localization  to  accelerate  through  the  remainder  of  the 
specimen. 

Vardoulakis  et  al  (1978)  discussed  the  spontaneous  formation  of  shear  bands  in 
sand  bodies.  The  formation  of  such  bands  is  looked  upon  as  a  bifurcation  problem.  They 
demonstrated  theoretically  and  experimentally  that  both  Coulomb’s  and  Roscoe’s  solutions 
for  the  orientation  of  the  shear  bands  are  possible,  namely, 

Oj.  =  jr/4  -  (p jl  (Coulomb) 

Sj.  =  ir/4  -  ^  Jl  (Roscoe) 
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where  ^  is  the  dilatancy  angle  (see  Arthur  and  Dunstan,  1982).  They  analyze  four  modes  of 
deformation  and  find  that  Coulomb’s  solution  holds  for  non-rotating  stress  axes  while  the 
Roscoe  solution  holds  for  co-rotating  stress  and  strain  increment  axes  during  failure. 

Vardoulakis  (1980)  in  a  continuation  of  the  previous  paper  gives  experimental 
results  obtained  in  biaxial  tests  on  sand  that  show  that  the  inclination  of  the  shear  bands 
indeed  fall  between  the  two  limits  of  Coulomb  and  Roscoe.  He  concludes  that  because  of 
the  failure  of  the  normality  rule  in  sands,  localization  under  those  conditions  always  occur 
in  the  hardening  regime. 

Vardoulakis  and  Graf  (1982)  placed  wooden  inclusions  in  their  biaxial  testing  device 
and  studied  shear  band  formation  as  an  imperfection  sensitivity  problem.  In  1983, 
Vardoulakis  studied  diffuse  and  localized  bifurcation  modes  in  axisymmetric,  rectilinear 
deformations  on  rigid-granular  dilatant  material.  He  found  that  localizations  in  the 
compression  test  occur  in  the  softening  regime  and  in  the  extension  test  they  occur  in  the 
hardening  regime. 

Hettler,  A.  and  Vardoulakis,  1.  (1984)  conducted  an  extensive  investigation  on  sands 
in  a  large  triaxial  apparatus  and  concluded  that  when  the  platens  were  lubricated 
bifurcation  into  diffuse  bulging  occurred  in  the  hardening  regime.  Diffuse  bulging  could  be 
avoided  by  using  flat  specimens.  However,  shear  bands  inevitably  occur  in  the  softening 
regime  of  the  stress-strain  curve. 

Vardoulakis  (1985)  examined  the  stability  and  bifurcation  of  water  saturated  sands 
and  normally  consolidated  clay.  He  described  their  behavior  by  means  of  a  two- 
dimensional  flow  theory  of  plasticity  with  a  non-associated  flow  rule.  He  concluded  that 
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for  dilatant  materials  the  dominant  failure  mode  is  shear  banding  that  occurs  close  to  the 
state  of  maximum  principal  effective  stress  ratio.  Whether  the  shear  bands  have  a 
thickness  or  not  is  a  topic  that  was  also  studied  by  Miihlhaus  and  Vardoulakis  (1987).  The 
conclusion  was  that  they  do  and  that  this  thickness  is  a  small  multiple  of  the  mean  grain 
size,  and  that  it  grows  as  one  progresses  in  the  strain  hardening  region. 

Having  proven  that  in  shear  band  analvsis  the  difference  between  the  true  stress 
rate  and  the  Jaumann  stress  rate  could  be  ignored,  Vermeer  (1982)  deduced  elegant  and 
simple  relations  predicting  the  instant  of  initiation  as  well  as  the  direction  of  the  shear 
bands.  His  predictions  are  in  line  with  Mandel’s  conclusions  for  material  stability. 
Molenkamp  (1985)  studied  the  material  behavior  under  uniform  deformation  and  the 
occurrence  of  a  shear  band  as  a  coupled  phenomenon.  He  attempted  to  solve  the 
bifurcation  problem  for  a  range  of  popular  constitutive  models  such  as  Mohr  Coulomb’s 
and  Lade  and  Duncan’s.  He  found  among  those  models  large  variations  in  the  predicted 
directions  and  instants  of  initiation  of  the  shear  bands. 

Desrues  and  Chambon  (1989)  examine  the  implication  of  using  complex  constitutive 
relations  on  the  development  of  shear  bands.  They  give  expression  of  the  restriction  to  be 
imposed  for  various  solutions  to  be  possible. 

Kolymbas  and  Rombach  (1989)  use  a  non-linear  constitutive  equation  for  sand  with 
fbced  materials  constants  to  express  analytically  a  necessary  criterion  for  shear  band 
formation.  This  band  is  considered  to  be  a  shearing  zone  which  may  appear  spontaneously 
in  the  course  of  a  homogeneous  plane  deformation.  Their  criterion  predicts  not  only  the 
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stress  state  at  the  moment  of  a  shear  band  formation  but  also  the  inclination  of  the  shear 


band  and  the  angle  of  initial  dilatancy  within  the  shear  zone. 

Bardet  (1991)  discusses  the  scatter  in  the  results  obtained  by  various  investigators 
while  emphasizing  the  observations  made  by  Desrues.  He  uses  an  extended  Mohr- 
Coulomb  model  to  get  a  more  accurate  analytical  expression  for  the  average  orientation  of 
the  shear  bands  in  soils.  The  concept  of  delayed  localization  together  with  this  model 
allows  him  to  retrieve  the  lower  and  upper  bounds  of  Roscoe  and  Mohr  Coulomb  for  the 
orientation  of  shear  bands. 

U  In  the  Areas  of  Strain  Localization  Using  Finite  Elements : 

There  are  two  main  approaches  to  treat  bifurcation  problems  in  the  context  of  finite 
elements  analysis;  the  direct  method  and  the  indirect  method.  The  direct  method  analyses 
the  bifurcation  condition  for  every  time  step,  for  each  element.  The  calculation  for  the 
bifurcation  branches  requires  the  determination  of  the  eigenvalues  and  eigenvectors  of  the 
tangent  stiffness  matrix.  In  the  indirect  method,  the  original  structure  is  slightly  modified 
with  the  help  of  initial  imperfections.  No  special  element  is  introduced  and  strain  softening 
is  often  part  of  the  model. 

Prevost  and  Hughes  (1981)  demonstrate  that  elastic-plastic  failure  states  may  be 
captured  in  finite  element  models  by  employing,  the  elastic-plastic  material  stiffiiess  to 
form  the  global  stiffness,  selective  integration  techniques  to  alleviate  mesh  locking  due  to 
incompressibility,  and  an  imperfection  in  the  form  of  weak  element.  Prevost  (1984) 
applied  this  technique  to  the  classical  punch  and  slope  stability  problems. 
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Bardet  and  Mortazavi  (1987),  introduced  a  material  imperfection  to  initiate  strain 
localization  in  a  homogeneously  strained  clay  subjected  to  plane  strain  compression. 

Hsu,  Peters  and  Saxena  (1987)  studied  the  influence  of  the  mesh  design  on  the 
location  of  initiation  of  shear  bands.  They  recommended  that  a  medium  size  mesh  be 
chosen  for  capturing  shear  bands. 

Potts,  Dounias  and  Vaughan  (1987)  used  the  direct  method  to  study  the  behavior  of 
frictional  materials  in  the  direct  shear  box.  The  boundary  conditions  themselves  initiate 
the  shear  bands.  Highly  stressed  zones  were  found  to  propagate  from  the  edges  of  the  box 
at  early  stages  of  loading.  Those  zones  grow  and  rotate  with  further  loading.  They 
examined  non-strain  softening  behavior  as  well  as  strain  softening  behavior  arising  from 
both  dilation  of  dense  soil  and  particle  orientation  in  clays.  In  this  last  case  the  analysis 
shows  that  at  the  peak  a  narrow  shear  zone  along  the  center  plane  of  the  sample  is  just 
beginning  to  form;  but  very  shortly  after  the  peak  the  shear  zone  is  completely  formed. 
The  post  peak  behavior  involves  plastic  straining  only  within  a  thin  central  line  of  elements 
with  elastic  unloading  of  the  material  on  either  side. 

Shuttle  and  Smith  (1988)  numerically  initiated  zones  of  localized  shear  strain  within 
idealized  cohesive  and  frictional  soil  specimens  by  assuming  inhomogeneities  and 
imperfect  boundaries.  For  non-strain  softening  cohesive  soil  the  localization  was  triggered 
by  material  imperfections.  For  strain  softening  cohesive  soils,  the  introduction  of  the 
softening  conditions  resulted  in  a  greater  or  lesser  degree  of  global  softening  or  brittleness. 
Upon  yield  major  stress  redistributions  take  place  as  localization  proceeds  through  the 
weakest  point  in  the  specimen.  Their  use  of  elastic  strain  until  the  peak  is  reached  results 
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in  the  storage  of  much  energy  which  is  later  used  to  drive  the  localization  process  on.  In 
reality  much  of  the  pre-peak  strain  in  soils  involves  dissipative  effects  and  a  smaller  amount 
of  energy  is  available  for  release.  The  authors  are  primarily  interested  in  pore  water 
pressure  migration  in  relation  to  shear  band  formation. 

Leroy  and  Ortiz  (1989)  confirm  that  in  the  absence  of  normality,  localization  can 
take  place  even  under  hardening  conditions.  Their  numerical  simulations  suggest  that  the 
observed  softening  is  a  geometrical  effect  not  attributable  to  constitutive  behavior  (see  also 
Drescher  and  Vardoulakis,  1982).  Conventional  finite  elements  techniques  are  shown  to  be 
limited.  They  propose  the  use  of  an  enriched  displacement  field  to  capture  the  localization 
at  the  point  of  bifurcation. 

Hobbs  and  Ord  (1989)  show  that  the  shear  band  width  is  sensitive  to  the  finite 
difference  mesh  chosen  in  the  solution  given  by  the  finite  difference  code  FLAG  and  is  also 
controlled  by  the  magnitude  of  both  angles  of  friction  and  dilation. 

1.4  In  the  Area  of  Fracture  Mechanics  Applied  to  Soils : 

Based  on  the  qualitative  descriptions  of  the  failure  mechanisms  responsible  for 
progressive  failure.  Palmer  and  Rice  (1973)  suggested  that  sliding  occurs  on  concentrated 
slip  surfaces.  Using  concepts  from  fracture  mechanics  (i.e.  the  J-integral),  attempts  were 
made  to  assess  the  time  dependence  governing  the  propagation  rate  of  a  particular  shear 
band.  The  assumed  model  used  for  the  shear  band  was  very  similar  to  the  cohesive  force 
models.  It  asserts  that  there  is  a  fixed  linear  relationship  between  the  shear  stress  in  the 
material  and  the  displacement  required  to  produce  it.  Since  displacement  is  the  integral  of 
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the  strains  over  a  given  region,  an  immediate  consequence  is  that  size  effects  will  occur. 
Specifically,  the  size  effect  appears  in  the  resistive  part  of  the  crack  driving  force  equation. 
In  addition.  Palmer  and  Rice  give  advice  for  experimentally  obtaining  the  energy  release 
rate  in  an  approximate  manner  by  computing  the  energy  imder  the  unloading  portion  of  the 
stress-strain  curve  in  a  particular  shear  box  test  Their  model,  however,  does  not  predict 
the  direction  of  propagation  of  the  shear  band  (crack  growth). 

Saada,  A.S.,  Chudnovsky,  A.  and  Kennedy,  M.R.  (1985)  applied  concepts  of  linear 
fracture  mechanics  to  overconsolidated  brittle  clays  tested  in  Modes  I  and  H.  They  used  a 
stability  criterion  to  predict  failure  of  infinite  slopes. 

Chudnovsky,  A.,  Saada,  A,S.  and  Lesser,  AJ.  (1988)  and  Lesser  (1989)  studied  the 
propagation  of  cracks  and  damage  zones  in  stiff  clays,  in  mode  H,  under  cyclic  loading. 
Both  macroscopic  and  microscopic  observations  led  to  the  conclusion  that  the  deformation 
in  the  material  primarily  occurs  at  the  level  of  the  clusters.  Such  clusters  coalesce  along 
shear  bands  on  both  sides  of  a  propagating  crack  and  form  the  damaged  zone.  Stress 
distributions  aroimd  the  tip  of  the  crack  were  also  given  by  Lesser. 

Using  a  plane  stress  direct  shear  apparatus  VaUejo  (1987)  tested  brittle  kaolinite 
with  pre-existing  opened  cracks.  He  found  that  the  cracks  propagated  from  the  tip  in  the 
direction  normal  to  the  direction  of  the  maximum  tensile  stress.  In  1988,  he  observed  the 
marked  differences  between  the  behavior  of  brittle  and  ductile  clay  when  subjecting 
specimens  with  inclined  cracks  to  compressive  loads.  For  the  particular  clay  used  (namely 
kaolinite)  he  found  that  at  a  water  content  less  than  20  percent  the  behavior  was  brittle. 
Vallejo  (1989)  also  examined  the  interference  between  two  cracks  in  a  specimen  of  brittle 
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clay  subjected  to  compression.  Depending  on  the  relative  position  of  the  cracks  and  the 
direction  of  the  load  he  noticed  the  weakening  effects  due  to  the  superposition  of  the 
tensile  stresses  that  develop  near  the  crack  tips. 
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CHAPTER  II 


EXPERIMENTAL  PROGRAM 


The  following  sections  describe  the  testing  and  recording  devices,  the  specimens 
preparation,  the  scope  of  the  investigation  and  the  conduct  of  the  tests. 

2.1  The  Hollow  Cylinder  Testing  Cell 

Used  since  1965  to  validate  failure  criteria  and  constitutive  equations,  this  cell  was 
first  used  in  1985  by  Saada  et  al.  to  study  fi-acture  and  damage  propagation  in  clay  soils.  It 
has  been  discussed  at  length  in  two  state-of-the-art  papers  by  Saada  and  Townsend  (1981) 
and  Saada  (1988).  Fig.  2.1  shows  the  cell  used  at  Case  Western  Reserve  University.  All 
the  measurements  are  made  both  inside  and  outside  the  cell  with  a  variety  of  transducers. 
Axial,  torsional  and  spherical  stresses  can  be  applied  to  the  specimen  through  a  piston  and 
a  confining  fluid.  To  maintain  uniformity  of  the  normal  stresses,  the  same  pressure,  namely 
the  cell  pressure,  acts  inside  and  outside  the  hollow  cylinder. 

Expressions  for  the  average  values  of  the  stresses  and  strains  in  terms  of  the 
measured  forces  and  displacements  are  given  in  the  two  papers  previously  mentioned. 

2.2  Loading  Devices 

Two  types  of  loading  devices  were  used  in  this  investigation: 

The  first  is  a  controlled  stress  device  which  applies  combinations  of  axial,  torsional 
and  spherical  stresses  at  predetermined  ratios,  thus  resulting  in  radial  paths.  It  is  called 
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SPAC  and  has  been  described  by  Saada  and  Townsend  (1981).  With  this  device  one 
cannot  go  beyond  the  peak,  and  failure  is  catastrophic  in  nature.  All  its  controls  are 
pneumatic  and  use  fluidic  components.  Fig.  2.2  shows  a  photograph  of  SPAC. 

The  second  is  a  controlled  rotation  device  that  uses  both  mechanical  and  fluidic 
controls.  A  motor  applies  a  slow  rotation  to  the  base  of  the  hollow  cylinder  cell  whose 
piston  is  prevented  from  rotating.  The  resulting  torque  is  measmed  with  a  transducer  and 
its  electrical  output  is  converted  to  a  pneumatic  one.  This  output  is  operated  on  by  a 
variety  of  fluidic  components  and  sent  to  an  actuator  which  applies  to  the  cell’s  piston  an 
axial  load  proportional  to  the  torque.  The  result  is  a  radial  path  with  the  possibility  of 
going  beyond  the  peak  when  the  failure  is  torsional  in  nature.  This  device  was  used 
whenever  one  wanted  to  observe  the  slow  development  of  the  shear  bands  and  the  slow 
propagation  of  damaged  zones  and  cracks.  Fig.  23  shows  a  photograph  of  this  device. 

All  the  tests  conducted  were  at  constant  volume.  If,  in  addition,  the  height  of  the 
specimen  is  fixed,  a  condition  of  constant  cross  section  is  present.  Fig.  2.4  shows  the  system 
of  stresses  and  the  sign  convention  adopted  in  this  paper. 

23  Observation  and  Recording  Devices 

Standard  and  high  speed  cameras  were  used  to  keep  the  visual  record  of  the  cracks 
and  shear  bands  propagation  in  the  test  specimens. 

A  Bioquant  surface  analyzer  was  used  to  study  the  fabric  of  the  material  on  a 
macroscopic  level.  This  device  helps  observe  and  quantify  changes  in  the  fabric,  specially 
in  the  vicinity  of  the  crack,  shear  bands  and  damaged  zones  (Fig.  2.5). 
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Transmission  electron  microscopes  were  used  to  look  inside  the  macroscopic  units 
examined  by  the  surface  analyzer. 

A  state-of-the-art  data  acquisition  and  reduction  system  was  used  to  obtain, 
calculate  and  graphically  present  the  data  (Fig.  2.6). 

2.4  Preparation  of  Clay  Specimens.  Macroscopic  Examination 

Two  kaolinitic  types  of  clay  were  used  in  this  investigation.  The  first  one  is  known 
as  Edgar  Plastic  Kaolin  (EPK).  It  has  been  extensively  used  in  soil  mechanics  research.  It 
has  a  low  activity,  a  liquid  limit  of  56.3%  and  a  plastic  limit  of  37.3%.  The  second  one  is 
known  as  Hydrite  121  (H121).  It  has  been  extensively  used  at  McGill  University,  Canada. 
It’s  liquid  limit  is  49.6%  and  its  plastic  limit  is  37.9%  resulting  in  a  low  plasticity  index.  It  is 
different  from  EPK  in  that  its  particle  are  quite  large  and,  when  dried  it  easily  crumbles. 
H 121  is  very  difficult  to  bring  to  a  low  water  content  by  consolidation  alone  and  one  has  to 
revert  to  slow  drying  to  obtain  stiff  specimens:  A  relatively  small  loss  of  water  causes  this 
clay  to  become  quite  stiff  compared  to  EPK. 

The  procedure  briefly  described  herein  has  been  presented  in  detail  by  Saada 
(1988).  The  clay  powder  was  mixed  under  vacuum  with  distilled  water  to  form  a  slurry  at 
approximately  twice  the  liquid  limit.  The  slurry  was  then  consolidated  in  a  large  20  cm 
diameter  consolidometer  to  yield  a  block  from  which  solid  cylinders  were  cut.  5cm 
diameter  cores  were  removed  from  the  center  of  the  cylinders  to  yield  the  specimen  shown 
schematically  in  Fig.  2.7a.  While  this  shape  was  most  satisfactory  for  routine  tests  not 
involving  crack  propagation  studies,  it  was  found  necessary  to  adopt  another  configuration 
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whenever  careful  of  observation  of  shear  bands  formations  were  needed.  The  specimens 
with  enlarged  ends  shown  in  Fig.  2.7b  were  made  by  pressing  and  coring  the  conlSguration 
shown  in  Fig.  2.7a  using  special  cutters  and  moulds.  End  effects  were  thus  considerably 
reduced,  if  not  totally  eliminated  during  testing.  Fig.  2.8  shows  a  succession  of  photographs 
illustrating  the  samples  preparation. 

When  studying  isotropic  clays,  the  powder  was  hand  mixed  and  specimens 
consolidated  isotropically. 

To  study  crack  propagation  radial  notches  were  made  in  the  specimens  using  a  thin 
vee-shaped  brass  sheet.  The  ratio  of  the  circumferential  length,  a,  of  the  notch  to  the  mean 
perimeter  1  of  the  specimen  varied  between  0.1  and  0.4.  Some  of  the  notches  were  inclined 
to  the  horizontal.  The  reference  ratio  was  a/1  =  0.18  (Fig.  2.7b). 

To  prevent  the  notch  (or  crack)  from  closing  during  consolidation  and  shear,  two 
layers  of  very  thin  Teflon  sheets  were  inserted  in  the  artificially  made  notch.  Such  sheets 
transmit  very  little  or  no  shear;  but  transmit  compressive  normal  stresses  across  the  notch 
surfaces. 

To  simplify  the  recording  of  the  observation  a  0.6  cm  square  grid  was  stamped  on 
the  specimens  (Fig.  2.9).  Photographs  were  taken  during  the  tests  to  study  the  deformation 
patterns.  Prior  to  placing  the  specimens  in  the  oven  for  moisture  content  determination, 
the  specimens  were  rotated  in  front  of  a  stationary  camera  and  photographs  were  taken  for 
small  increments  of  rotation.  Those  photographs  were  then  overlapped  with  matching 
grids  so  that  deformations  and  slip  patterns  could  be  examined  on  the  development  of  the 
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cylinder.  In  case  of  formation  of  shear  bands,  a  digitizer  was  used  to  obtain  a  line  drawing 
of  the  bands. 

Most  of  the  tests  were  conducted  on  consolidated  specimens.  This  kind  of 
consolidation  is  routinely  performed  using  a  pneumatic  analog  computer  that  measures  the 
water  expelled  from  the  specimen  as  the  pressure  in  the  cell  increases  and  imposes  a 
vertical  displacement  such  that  the  original  cross  section  remains  constant.  The  value  of 
Kq  is  continuously  given  by  a  transducer  during  the  whole  process.  At  the  end  of  the  Kq 
consolidation  process,  Kq  was  0.44  for  the  EPK  clay  and  0.51  for  the  H121  clay. 

All  the  specimens  were  rebounded  to  a  condition  of  spherical  state  of  stress  prior  to 
shearing.  Any  consolidation  test  in  which  the  piston  load  is  removed,  results  in  a 
specimen  that  is  normally  consolidated  in  the  horizontal  direction  and  that  has  an 
overconsolidative  ratio  (OCR)  of  nearly  2  in  the  vertical  direction.  If  one  is  to  think  in 
terms  of  mean  stresses  this  ratio  is  1.3.  Because  of  the  size  of  this  ratio  the  material  will  be 
referred  to  as  slightly  overconsolidated  or  normally  consolidated.  If  in  addition  to 
removing  the  piston  load  one  decreases  the  pressure  in  the  cell  by  a  factor  of  3  say,  the 
specimen  is  at  an  OCR  of  3  laterally  and  6  vertically.  In  terms  of  mean  stresses  the  OCR  is 
4.  The  return  to  a  spherical  state  of  stress  allows  one  to  conduct  a  radial  path  test. 

To  examine  the  clay  in  the  surface  analyzer  thin  sections  of  wax  impregnated 
specimens  are  needed.  At  the  end  of  a  test,  triaxial  specimens  are  placed  in  an  oven  to 
determine  the  final  water  content.  Once  this  is  done,  a  piece  of  the  dried  specimen 
including  the  zone  to  be  examined  is  placed  in  a  bath  of  wax  under  a  vacuum  bell.  It  is 
kept  there  for  a  specific  length  of  time  then  withdrawn  and  left  to  dry.  Thin  radial  sections. 
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3mm  thick  are  cut,  their  faces  polished,  and  then  placed  in  the  analyzer’s  microscope.  The 
image  that  appears  on  the  screen  lends  itself  to  a  variety  of  me^urements  controlled  by  the 
Bioquant  software. 

It  is  extremely  time  consiuning  to  find  an  impregnation  technique  that  yields  good 
resolution,  good  thin  sections  and  a  satisfactory  visual  demonstration  of  the  changes  in  the 
fabric  of  the  clay.  Through  trial  and  error  a  technique  to  impregnate  dry  specimens  of  clay 
with  wax  under  vacuum  was  developed.  Fig.  2.11  shows  a  section  in  an  undisturbed  part  of 
a  clay  specimen.  The  white  clusters  of  clay  are  surrounded  by  a  dark  matrix  of  wax  and  are 
homogeneously  distributed.  The  impregnation  technique  was  adjusted  such  that  clay 
clusters  of  about  100  microns  are  formed  and  are  visible  with  relatively  small 
magnification. 

.  Electron  microscope  examinations  were  made  at  McGill  University,  Canada,  by  Dr. 
R.  Yong. 

2.5  Scope  of  the  Investigation 

This  experimental  investigation  is  aimed  at  studying  the  mechanism  of  deformation 
and  the  development  of  the  shear  bands  in  a  clay  subjected  to  various  systems  of  stresses. 
Emphasis  is  placed  on  the  influence  of  cracks,  the  way  they  propagate  together  with  the 
shear  bands,  and  the  changes  that  occur  in  the  fabric  of  the  material,  in  particular  in  the 
vicinity  of  the  cracks  and  the  damaged  zone. 

Series  of  tests  were  conducted  on  specimens  with  the  configurations  in  Figs.  2.7a 
and  2.7b.  Most  of  the  cracks  artificially  induced  in  the  specimens  were  horizontal.  Only  a 
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few  were  inclined.  Some  specimens  had  two  cracks.  The  specimens  prepared  through  one 
dimensional  consolidation  are  cross  anisotropic  in  nature.  Some  tests  were  conducted  on 
isotropic  materials  to  examine  the  influence  of  anisotropy  on  the  direction  of  the  shear 
bands.  A  number  of  cyclic  loading  and  fast  loading  tests  were  also  conducted.  All  the  tests 
were  of  the  consolidated  undrained  type  with  pore  pressures  measurements. 

In  all,  the  experimental  investigation  was  quite  extensive.  Reference  tests,  i.e.,  tests 
without  cracks  had  to  be  conducted  to  allow  for  comparisons.  Their  number  and  variety 
are  sufficient  to  provide  a  data  base  to  validate  any  constitutive  model  for  cohesive  soils. 
The  tests  conducted  on  specimens  with  cracks  can  be  used,  not  only  to  describe  the 
kinematics  of  the  deformation,  but  also  to  test  the  appropriateness  of  presently  available 
crack  and  bifurcation  theories.  Sections  made  in  the  various  specimens  allow  one  to 
observe  how  the  fabric  changes. 

A  large  number  of  backup  tests  were  conducted  to  insure  repeatability.  In  addition 
one  dimensional  as  well  as  triaxial  consolidation  tests  were  needed  to  determine  properties 
and  provide  a  complete  description  of  the  clays  used.  Resonant  column  tests  on  specimens 
with  and  without  cracks  were  conducted  to  study  the  influence  of  the  cracks  on  shear 
moduli  and  damping  ratios.  All  the  tests  will  be  referred  to  in  the  following  discussions  and 
are  completely  detailed  in  the  appendices. 

Tables  2.1  to  2.9  list  the  tests  conducted  and  have  comments  giving  information  on 
the  specific  tests.  The  column  labelled  DESIGNATION  is  a  combination  of  a  number  and 
a  date.  Such  labels  are  used  to  locate  the  file  containing  the  data.  The  ang\3  is  the 
inclination  of  the  principal  stresses  on  the  vertical  axis  of  symmetry.  In  all  the  tests  for 
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which  p  is  given,  the  stress  path  was  radial.  To  each  /9  there  corresponds  a  value  for  the 
coefficient  b.  The  tests  m  which  the  length  is  kept  constant,  have  a  variable  ratio  of  axial  to 
torsional  stresses.  The  consolidation  pressure  in  the  tables  is  expressed  in  terms  of 
effective  cell  pressure  a,.;  with  the  understanding  that  the  specimen  has  been  Kq 
consolidated  at  a  cell  pressure  and  at  whatever  vertical  stress  is  necessary  to  keep  the 
cross  section  constant.  A  sketch  shows  the  configuration  used. 

The  order  of  the  tables  is  such  the  tables  listing  the  static  tests  on  EPK  are  given 
first,  followed  by  those  listing  the  dynamic  tests.  Also,  those  related  to  overconsolidated 
specimens  follow  those  related  to  Kg  normally  consolidated  specimens.  The  tests 
conducted  on  H121  clay  are  in  the  last  table. 

Table  2.1  lists  three  series  of  five  tests  on  EPK,  each  series  corresponding  to  a 
different  consolidation  cell  pressure  a^. ,  namely  241.5  kPa  (35  psi),  345  kPa  (50  psi)  and 
552  kPa  (80  psi). 

Table  2.2  lists  three  tests  on  EPK  used  to  examine  the  shear  bands  that  develop  in 
specimens  without  cracks.  Tv/o  tests  were  conducted  under  controlled  rotation  conditions 
and  one  under  controlled  stress  condition. 

Table  2.3  lists  test  conducted  on  EPK  on  specimens  with  horizontal  cracks  xmder 
various  inclinations  of  principal  stresses,  five  free  to  change  It  th  and  two  with  fixed 
length;  also  listed  are  three  tests  with  inclined  cracks. 

Table  2.4  lists  tests  conducted  on  K^  overconsolidated  EPK  where  the  cell  pressure 
Uj.  dropped  from  621  kPa  (90  psi)  to  207  kPa  (30  psi)  with  corresponding  drop  in  the 
vertical  effective  stress  as  explained  in  the  previous  section.  There  were  no  initial  cracks  in 
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those  spedmens.  They  are  used  as  reference,  as  well  as  to  study  shear  bands  development 
with  the  inclination  of  the  principal  stresses. 

Table  2.5  lists  tests  similar  to  the  ones  appearing  in  Table  4  but  with  a  horizontal 
crack  whose  ratio  a/1  is  equal  to  0.18.  In  addition,  it  lists  tests  conducted  in  torsion  on 
specimens  with  horizontal  cracks  with  different  ratios  a/1. 

Table  2.6  lists  two  torsional  cyclic  tests  and  one  torsional  impulse  test  conducted  on 
EPK  with  Kq  consolidated  specimens  (5^  =  345  kPa)  and  with  fixed  length.  One  of  the 
specimens  tested  under  <yclic  condition  had  a  45°  inclined  crack  and  the  two  others  had 
horizontal  cracks. 

Table  2.7  lists  two  torsional  cyclic  tests  conducted  on  EPK  with  Kq  overconsolidated 
specimens,  both  with  horizontal  cracks,  one  with  fixed  length  and  the  other  free  to  change 
length. 

Table  2.8  lists  five  tests  conducted  on  EPK  with  K^  overconsolidated  specimens  and 
two  cracks.  The  vertical  distance  between  the  cracks  d  is  0.63  cm  (0.25  in)  and  the 
horizontal  distance  e  varies  between  zero  and  2.54  cm  (1.0  in.).  Four  of  the  specimens 
have  their  cracks  in  the  downstep  relative  position  and  one  in  the  upstep  relative  position. 
All  were  free  to  deform  longitudinally. 

Table  2.9  lists  all  the  tests  conducted  on  the  H121  clay.  There  are  five  reference 
static  tests  with  various  inclinations  of  principal  stresses,  two  torsion  tests,  and  one  fixed 
length  torsional  cyclic  tests. 

Not  listed  in  the  tables  are  the  limited  number  of  the  resonant  column  tests  and  the 
identification  triaxial  and  consolidation  tests. 
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TABLE  2.1  -  STATIC  TESTS  ON  Kq  CONSOLIDATED  EPK 
CONFIGURATION  OF  HG.  2.7a,  NO  CRACK 


po 

— 

b 

a(.kPa 

_ 345 _ 

DESIGNATION 

0 

0 

31-1-29-91 

lb-1-10-90 

36-2-15-91 

30 

.25 

32-1-31-91 

2-11-15-89 

37-2-19-91 

45 

.50 

33-2-5-91 

3b- 12-28-89 

38-2-21-91 

60 

.75 

35-2-12-91 

5b- 11-13-91 

40-3-5-91 

90 

1.0 

34-2-7-91 

4b- 12-20-89 

41-3-7-91 

COMMENTS  AND 
CONFIGURATION 


.Stress  Controlled 
.Radial  Stress  Path 
.Free  Axial  Deformation 


Average  w% 


36.20 


34.53 


32.66 


TABLE  2.2  -  STATIC  TESTS  ON  Kq  CONSOLIDATED  EPK 
CONFIGURATION  OF  FIG.  2.7b,  NO  CRACK 


- -  1 

COMMENTS  AND 

po 

b 

Of.  kPa 

DESIGNATION 

CONTROL 

CONFIGURATION 

30 

.25 

■i 

62-5-16-91 

ROTATION 

.Free  Axial  Deformation 

45 

.50 

Ha 

3t-l-23-90 

STRESS 

60 

.75 

345 

61-5-14-91 

ROTATION 

-1 

PT - 

a 

- 

60 

.75 

345 

61b-19-ll-91 

STRESS 

1 

Average  w%  =  34.42 

— 

L_ 

5 

- 
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TABLE  23  -  STATIC  TESTS  ON  CONSOLIDATED  EPK,  =  345  kPa 


CONHGLnRATION  OF  FIG  2.7b,  WITH  CRACK 


po 

b 

DESIGNATION 

CONTROL 

CRACK 

ANGLE 

COMMENTS  AND 
CONHOURATION 

0 

30 

45 

60 

90 

0 

.25 

.50 

.75 

1.0 

504-3-91 

514-5-91 

6V-9-26-90 

524-9-91 

534-11-91 

STRESS 

STRESS 

ROTATION 

STRESS 

STRESS 

0 

0 

0 

0 

0 

.Free  Axis 
.a/1  =  .18 

dDef 

ormation 

TORSION 

3f-2-6-90 

6f-6-2190 

STRESS 

ROTATION 

0 

0 

.Fixed  Length 
.a/U.18 

1 

564-17-91 

574-19-91 

584-23-91 

ROTATION 

ROTATION 

ROTATION 

45»/ 
45“  \ 
22.5 

.Free  Axial  Deformation 

1 

/ 

K 

1 

il 

i 

Average  w%  =  34.60 

— 

1 

>  1- 
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TABLE  2.4  -  STATIC  TESTS  ON  Ko  OVERCONSOLIDATED  EPK 
MEAN  OCR  =  4,  CONFIGURATION  OF  FIG  2.7b,  NO  CRACK 


po 

b 

DESIGNATION 

CONTROL 

0 

0 

87-7-25-91 

STRESS 

30 

.25 

89-8-1-91 

ROTATION 

45 

.50 

86-7-23-91 

ROTATION 

60 

.75 

90-8-13-91 

ROTATION 

90 

1.0 

88-7-30-91 

STRESS 

COMMENTS  AND 
CONHGURATION 


.Free  Axial  Deformation 


Average  w%  =  33.92 


•  One  static  test  on  EPK,  DESIGNATION  7a-9-6-90  was  conducted  on  an  isotropic 
overconsolidated  specimen  under  controlled  deformation  conditions,  and  fixed  len^ 
under  torsion.  OCR  =  3,  w%  =  30.20. 
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TABLE  2^  -  STATIC  TESTS  ON  Ko  OVERCONSOLIDATED  EPK 
MEAN  OCR  =  4,  CONHOURATION  OF  FIG.  2.7b,  WITH  CRACK 


po 

B 

DESIGNATION 

CONTROL 

a/1 

COMMENTS  AND 
CONFIGURATION 

0 

30 

45 

60 

90 

_ 

0 

.25 

.50 

.75 

1.0 

67- 5-31-91 

68- 6-5-91 
66-5-29-91 
70-6-14-91 

69- 6-7-91 

STRESS 

ROTATION 

ROTATION 

ROTATION 

STRESS 

.18 

.18 

.18 

.18 

.18 

.Free 

_i 

Axial  D( 

! 

efonnation 

L 

TORSION 

19-1-23-91 

16-12-13-90 

18-1-15-91 

15-11-28-90 

13- 11-7-90 

14- 11-20-90 

ROTATION 

ROTATION 

ROTATION 

ROTATION 

ROTATION 

ROTATION 

.Fixec 

■7 

/ 

/ 

/ 

\  Length 

! 

Average  w%  =  32.97 
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TABLE  2.6  -  DYNAMIC  TESTS  ON  CONSOLIDATED  EPK,  =  345  kPa 
CONHOURATION  OF  FIG.  2.7b,  WITH  CRACK 


CRACK 

b 

DESIGNATION 

CONTROL 

ANGLE 

TORSION 

3-1-30-90 

CYCUC 

0 

3-2-13-90 

IMPULSE 

0 

3-2-22-90 

CYCUC 

45 

COMMENTS  AND 
CONFIGURATION 


.Fixed  Length 
.a/1  =  .18 


Average  w%  =  33.87 


TABLE  2.7  -  DYNAMIC  TESTS  ON  Kq  OVERCONSOUDATED  EPK 
MEAN  OCR  =  4,  CONHGURATION  OF  FIG.  2.7b,  WITH  CRACK 


0 

b 

DESIGNATION 

CONTROL 

a/1 

COMMENTS  AND 
CONHGURATION 

TORSION 
45  1  .5 

17-1-3-91 

6VC-10-3-901 

CYCUC 

CYCUC 

.1 

.18 

.Fixed 
.Free  i 
T/ 

I 

th 

Deformation 

f 

Average  w%  =  32.03 

•  One  cyclic  loading  test  on  EPK,  DESIGNATION  7-6-27-90  was  conducted  on  an 
isotropic  overconsolidated  specimen  with  fixed  length  under  torsion.  OCR  =  3,  w%  = 
29.91% 
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TABLE  2.8  -  STATIC  TESTS  ON  Ko  OVERCONSOUDATED  EPK 
MEAN  OCR  =  4,  WITH  TWO  CRACKS,  FREE  AXIAL  DEFORMATION 


p 

■ 

DESIGNATION 

CONTROL 

CRACKS 

POSITION 

COMMENTS 

45 

.5 

81-7-3-91 

ROTATION 

! _ j_ 

i  d 

.Downstep 
e  =  0  d  =  .6cm 

45 

.5 

82-7-9-91 

ROTATION 

“  T 

1 

i 

a/1  =  .18 

.Upstep 

e  =  0  d  =  .6cm 

45 

.5 

83-7-11-91 

ROTATION 

a/1  =  .18 

.Downstep 
e  =  .6cm  d  =  .6cm 

45 

.5 

84-7-17-91 

ROTATION 

1 

1 

-jeK 

1 :  ,1— 

a/1  =  .18 

.Downstep 
e  =  1.3cm  d  =  .6cm 

45 

.5 

85-7-19-91 

ROTATION 

a/1  =  .18 

.Downstep 
e  =  2.54cm  d=  .6cm 

1  !  1 

a/1  =  .18 

Average  w%  =  33.62 


TABLE  2.9  -  TESTS  ON  Kq  CONSOLIDATED  H121,  =  345  kPa 


TORSION 

TORSION 


TORSION 


DESIGNATION 

CONTROL 

1-5-1-90 

STRESS 

2-6-1-90 

STRESS 

3-5-3-90 

STRESS 

5-6-8-90 

STRESS 

4-6-5-90 

STRESS 

3f-5-16-90 

STRESS 

6f-9-17-90 

ROTATION 

3-6-13-90 

CYCUC 

COMMENTS  AND 
CONFIGURATION 


.Free  Axial  Deformation 


.Fixed  Length 
.Crack  a/1  =  .18 


Average  w%  =  31.14 


MAROCNCD  S,S. 


BRASS  SLEEVE 


NOOIZEO  AL. 


0  RING 

transducer 


DRAINA6E 


y-— RLEXI6LASS 

S.S.  TOOTHED 
POROUS  STONE 

SPECIMEN 


1  for  Axial  and  Torsional  Stresses  on  Hollow  Cylinders 
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Figure  12  Stress  Controlled  SPAC  and  Supporting  Equipment 


Figure  23  Defomiatioii  Controlled  SPAC  and  Supporting  Equipment 
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Figure  2.4  Notation 
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automatic  data  acquisition  system 
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Figure  2.6  Schematic  Set-up  for  Hollow  Cylinder  System 


Figure  2.8  Preparation  of  a  Hollow  Cylinder  of  Clay 
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Figure  2.9  Specimens  with  Square  Grid 


VACUUM 


Schematic  for  Specimen  Impregnation. 


Figure  2.10  Impregnation  and  Sectioning  of  Clay  Specimens 
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UNDAMAGED 

REGION 


Figure  2.11  Clay  Clusters  after  Impregnation 
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TEST  RESULTS  AND  OBSERVATIONS  RELAIED  TO  STRENGTH 

3.1  Observations  Corresponding  to  the  Tests  in  Tables  2.1, 22  and  2S, 

The  test  conducted  on  specimens  without  cracks  serve  not  only  as  a  reference  in  the 
analysis  of  strength  and  shear  band  directions,  but  also  to  determine  the  influence  of  the 
cracks  on  the  strength  parameters  as  the  system  of  stresses  rotates. 

Figs.  3.1a  to  3.1e  show  the  normalized  octahedral  shearing  stresses  versus  the 
octahedral  shearing  strains  various  consolidation  pressures  and  various  inclinations 

^  01  khe  principal  stresses.  Notice  that  b  =  (02  •  03)7(01  -  03)  =  sin2/9.  The  octahedral 
variables  take  into  account  the  influence  of  03  and  by  comparing  the  various  figures  one 
can  see  the  influence  of  anisotropy.  Regardless  of  the  consolidation  pressure  the  material 
behaves  in  a  stiffer  way  when  the  compressive  component  of  the  applied  stress  is  prevalent. 
This  is  primarily  due  to  the  transverse  isotropy  resulting  from  Kj,  consolidation.  As  the 
strain  increases  the  slope  of  the  stress-strain  curves  decreases  far  more  rapidly  when  the 
tensile  stresses  are  the  predominant  ones.  For  the  same  direction  of  the  principal  stresses 
normalization  does  not  lead  to  a  unique  curve  as  can  be  seen  in  each  of  the  figures.  The 
situation  is  made  more  dramatic  when  one  tries  to  superimpose  all  the  stress  strain  curves. 
Uniqueness  is  certainly  not  present.  For  completeness  the  pore  water  pressures  are  given 
in  Figs.  3.2a  to  3.2e.  Notice  that  in  those  figures  the  value  of  the  applied  mean  stress  is 
subtracted  from  the  developed  pore  pressure  and  the  result  normalized  with  respect  to  the 
effective  consolidation  pressure.  Here,  too,  this  normalized  pore  water  pressure  due  to 
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shear  varies  dramatically  with  the  inclination  of  the  principal  stresses.  The  above  puts  in 
question  the  theories  based  on  unique  pore  water  pressure  coefficients  at  failure,  be  they 
Skempton’s  or  Henkel’s. 

Figs.  3.3a  and  3.3b  gives  the  results  obtained  for  the  H121  clay.  Here,  again, 
uniqueness  is  not  obtained. 

To  assess  the  validity  of  various  strength  theories,  it  is  interesting  to  study  the 
direction  of  the  shear  bands  that  were  observed  in  the  specimens  when  subjected  to  various 
stress  paths.  There  are  five  different  stress  paths  and,  if  the  Mohr  Coulomb  theory  of 
failure  was  perfectly  valid,  the  plane  of  failure  at  each  point  should  contain  the 
intermediate  principal  stress,  which,  in  our  case  is  radial.  The  tests  with  paths 
corresponding  to  an  inclination  0  =  30, 45  and  60  degrees  show  that  this  is  not  the  case: 

The  surfaces  in  the  three  cases  are  inclined  on  the  radial  direction,  the  most  noticeable 
being  the  ones  corresponding  the  0  =  30°  tests.  In  this  case,  however,  the  inclination  was 
under  ten  degrees.  Because  of  symmetry  the  pure  compression  and  extension  tests  give  the 
standard  conical  failure  surface.  Fig.  3.4  shows  the  angle  of  friction  <p '  as  defined  by  sin'^ 
(ci  -  03)7(01  +  03)  obtained  fi’om  the  tests  conducted  on  EPK  as  a  function  of  the 
inclination  0  (or  b).  We  see  that  ^  *  varies  drastically  with  0  and  loses  its  meaning  as  the 
tensile  stresses  become  predominant.  The  variations  in  <p*  among  the  various  one¬ 
dimensional  consolidation  pressures  are  primarily  due  to  changes  in  anisotropy  as  the 
consolidation  pressures  are  increased.  Similar  results  are  obtained  for  H121  illustrating  the 
influence  of  both  b  and  anisotropy  on  the  friction  angle  0'.  Such  observations  have 
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appeared  in  the  literature  and  commented  on  by  Saada  and  Bianchini  (1975)  and  Ladd  et 
al.  (1977). 

It  is  worthwhile  to  comment  on  the  tests  listed  in  Table  22.  These  tests  were 
conducted  on  the  same  material  used  in  Table  2.1  but  with  a  different  geometric 
configuration  and  in  two  cases  at  a  controlled  rate  of  rotation  rather  than  at  a  controlled 
rate  of  stress.  Figs.  3.5  shows  the  superposition  of  the  results  on  the  ones  obtained  from 
the  tests  listed  in  Table  2.1.  The  corresponding  curves  are  quite  close,  the  difference  being 
within  acceptable  experimental  error.  In  addition,  it  is  appropriate  to  show  at  this  stage  the 
inclination  of  the  shear  bands  for  those  three  tests.  As  previously  stated,  the  photographs 
in  Figs.  3.6a  and  b  were  made  of  a  succession  of  single  shots  taken  while  rotating  the 
specimen  around  its  axis  of  symmetry  prior  to  placing  it  in  the  oven  for  moisture  content 
determination.  Many  photographs  were  also  taken  during  the  test  to  observe  the 
development  of  the  shear  bands.  One  notices  on  those  photographs  that  the  inclination  of 
the  shear  bands  varies  as  one  moves  from  compression  to  extension.  Shown  on  the 
digitized  plots  3.6az  and  3.6bz  is  the  direction  of  shear  bands  predicted  by  the  Mohr- 
Coulomb  theory.  The  angle  (45  +  0 '  /2)  giving  the  inclination  of  the  Coulomb  slip  planes 
was  computed  with  an  average  angle  of  friction  of  33  degrees  obtained  from  direct  triaxial 
compression  tests  (Fig.  3.4).  One  can  see  that  for  the  two  tests  at  hand  the  agreement  is 
remarkable.  The  grid  stamped  on  the  specimens  allows  us  to  measure  angles  with  quite  a 
reasonable  accuracy.  It  also  visually  enhances  the  shear  bands  and  their  directions. 

While  the  observations  above  are  important  and  are  the  first  of  their  kind  to  be 
obtained  under  combined  stresses,  the  fact,  remains  that  the  reference  <p '  was  obtained 
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from  direct  compression  and  Fig.  3.4  is  enough  of  a  reminder  that  the  traditionally  defined 
^ '  is  certainly  not  unique. 

3.2  Observations  Corresponding  to  the  Tests  in  Table  23. 

Seven  out  of  the  ten  tests  listed  in  Table  23  were  conducted  on  specimens  with 
horizontal  cracks  and  a  ratio  a/1  =  0.18.  The  three  last  tests  had  inclined  cracks  2.5  cm  in 
length.  Figs.  3.7a  and  b  show  the  stress-strain  curves  for  various  inclinations  of  the 
principal  stresses  on  the  axis  of  symmetry  and  the  pore  water  pressures  respectively.  In 
Figs.  3.7  we  have  superimposed  the  curves  from  Fig.  3.5  to  study  the  influence  of  the  crack. 
Such  influence  is  seen  to  be  minimal  for  the  size  of  the  crack  studied.  It  will  be  seen  in  the 
chapter  reserved  to  the  kinematics  of  deformation  that  once  the  shear  band  starts  in  the 
vicinity  of  a  crack  the  deformation  is  completely  controlled  by  the  plastic  behavior  of  the 
material. 

In  Figs.  3.8a  and  b  the  results  of  three  torsion  tests  on  specimens  with  horizontal 
cracks  are  superimposed;  two  conducted  under  fixed  length  conditions  and  one  free  to 
deform  axially.  In  the  fixed  length  cases  the  principal  stresses  rotate  a  little  depending  on 
the  axial  stress  that  is  developed  during  torsion.  The  presence  of  the  axial  stress  necessary 
to  prevent  elongation  results  in  higher  stresses  for  the  same  strain  and  in  differences  in  the 
measured  pore  water  pressures. 

The  influence  of  the  inclination  of  the  crack  on  the  shear  strength  behavior  and  the 
pore  water  pressure  is  shown  in  Figs.  3.9a  and  b.  There,  all  the  specimens  are  under  a  state 
of  pure  torsion  and  free  to  deform  axially.  Whether  the  inclination  is  45,  22.5  or  0  degrees 
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the  undrained  strength  is  the  same;  except  at  veiy  large  strains.  The  pore  water  pressures 
are  close  to  each  others.  Here  again,  it  appears  that  once  the  shear  bands  develop  the 
behavior  is  plastic  and  the  influence  of  the  crack  nearly  vanishes.  This  will  be  further 
discussed  in  the  next  chapter. 

33  Observations  Corresponding  to  the  Tests  in  Table  2.4. 

This  set  of  tests  was  conducted  on  overconsolidated  specimens  of  EPK  clay  without 
cracks  and  a  mean  overconsolidation  ratio  of  4.  Like  their  counterparts  in  Table  2.1,  they 
serve  as  a  reference;  this  time  for  the  study  of  strength,  crack  and  damage  propagation  in  a 
brittle  material.  Figs  3.10a  and  b  show  the  strength  and  pore  water  pressure  relations.  Fig. 
3.11  a,  b  and  c  shows  the  pattern  of  slip  lines  that  developed  for  an  angle  p  =  30°,  )8  =  45° 
and  j8  =  60°.  In  all  cases  the  prediction  of  the  inclinations  of  the  slip  lines  by  Coulomb’s 
Theory  (with  0 '  equal  to  33  degrees)  fits  extremely  well  the  measured  ones;  as  shown  on 
the  digitized  plots  3.11az,  3.11  bz  and  3.11cz.  The  photographs  in  those  figures  will  be  used 
for  comparison  purposes  when  examining  the  shear  bands  that  develop  in  the  presence  of  a 
crack. 

3.4  Observations  Corresponding  to  the  Tests  in  Table  2.5. 

There  are  two  sets  of  tests  in  this  table:  One  set  in  which  the  specimens  with  an  a/1 
=  0.18  were  subjected  to  the  same  radial  stress  paths  listed  in  Table  2.4,  and  a  second  set  in 
w...:h  the  specimens,  with  various  sizes  of  cracks,  were  subjected  to  torsion  under  fixed 
length  condition.  Figs.  3.12a  and  b  show  the  shear  strength  and  the  pore  water  pressure  as 
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a  function  of  the  strains.  In  these  figures  we  have  superimposed  the  results  shown  in  Figs. 
3.10  to  study  the  influence  of  the  crack.  As  in  the  case  of  Figs.  3.7  such  an  influence  is  seen 
to  be  minimal  in  spite  of  the  length  of  the  crack  being  about  twenty  percent  of  the 
perimeter. 

The  results  above  led  to  a  study  of  the  influence  of  the  length  of  the  crack  on  the 
strength.  Six  specimens  with  various  values  of  a/1  were  subjected  to  torsion  under  fixed 
length  conditions.  The  relations  stress  versus  strain  and  pore  water  pressure  versus  strain 
are  shown  in  Figs.  3.13a  and  b.  The  ultimate  strength  drops  as  the  size  of  the  crack 
increases.  This  is  expected.  However  the  fact  that  the  length  remains  constant  results  in 
higher  values  of  the  stresses  when  comparisons  are  made  with  tests  in  which  the  specimens 
are  fi’ee  to  elongate. 

The  photographs  taken  during  the  formation  of  the  shear  bands  in  the  vicinity  of  the 
tip  of  the  crack  will  be  examined  in  Chapter  IV  in  a  comparative  study  of  the  kinematics  of 
specimens  with  and  without  cracks. 

3.5  Observations  Corresponding  to  the  Dynamic  Tests  in  Tables  2.6  and  2.7. 

The  tests  in  those  tables  were  conducted  primarily  to  examine  the  pattern  of 
damage  propagation  when  a  specimen  with  a  crack  is  subjected  to  cyclic  and  impulse 
loadings.  Photographs  were  taken  as  torques  of  fixed  amplitude  were  cyclically  applied  to 
the  specimen,  and  as  a  step  loading  was  induced.  The  third  test  in  Table  2.6  was  conducted 
on  a  specimen  with  an  inclined  crack.  Table  2.7  lists  the  dynamic  tests  conducted  on 
overconsolidated  material. 
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Fig.  3.14  shows  the  decay  of  the  shear  modulus  and  increase  in  damping  for  the  first 
specimen  listed  Table  2.6.  The  hyterisis  loops  increase  in  size  as  their  inclination  on  the 
horizontal  axis  decreases  with  the  number  of  <ycles.  Note  the  lack  of  symmetry  in  Fig.  3.14. 
On  the  other  hand  Fig.  3.15  which  corresponds  to  the  isotropically  overconsolidated 
specimen  shows  a  more  symmetric  pattern.  Such  hysterisis  loops  had  been  noticed  and 
commented  on  by  Macky  and  Saada  (1984)  and  Saada  and  Macky  (1985)  for  both 
laboratory  prepared  specimens  and  undisturbed  specimens. 

Fig.  3.16  indicates  that  the  hysterisis  loops  for  an  inclined  crack  in  a  Kq  consolidated 
material  are  not  different  from  the  ones  obtained  with  a  horizontal  one  for  an  isotropic 
clay.  This  shows  that  the  fabric,  the  inclination  of  the  crack  and  the  direction  of  the  stress, 
all  affect  the  response. 

The  two  tests  in  Table  2.7  yield  results  similar  to  the  ones  above.  The  development 
of  shear  bands  in  all  the  tests  listed  in  Tables  2.6  and  2.7  was  carefully  observed  and 
photographed.  Crack  and  damage  propagations  were  not  expected  in  the  K^,  normally 
consolidated  specimens;  but  they  did  take  place,  just  as  they  did  with  overconsolidated 
specimens.  Comparisons  will  be  made  in  Chapter  IV. 

3.6  Observations  Corresponding  to  the  Tests  with  Two  Cracks  in  Table  2.8. 

Those  tests  were  conducted  to  answer  questions  about  coalescence  of  cracks  and 
their  damage  zone  during  the  shearing  process.  In  the  five  tests  listed  in  Table  2.8  the  two 
cracks  were  separated  horizontally  and  vertically  as  shown.  While  d  was  kept  constant  at 
0.6  cm  (0.25  in),  e  varied  between  0  and  2.54  cm  (1.0  in.).  The  direction  of  the  torque  was 
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the  same  for  all  tests  and  the  specimens  free  to  deform  axially  leading  to  the  same  value  of 
/8  =  45°.  Notice  that  the  two  first  tests  are  identical  except  that  the  first  has  its  cracks  in 
the  downstep  position  while  the  second  has  them  in  the  upstep  position. 

Figs.  3.17a  and  b  shows  the  strength  and  pore  pressure  curves  for  all  five  tests  in 
Table  2.8.  Also,  drawn  on  this  figure  for  comparison  purposes  is  the  curves  for  a  specimen 
with  one  crack  from  Fig.  3.12,  and  the  curve  for  a  specimen  with  no  cracks  from  Fig.  3.10. 
As  expected  the  strength  curves  for  two  cracks  fall  below  the  ones  with  one  or  no  crack. 
Also,  the  position  of  the  cracks  relative  to  the  direction  of  the  shearing  stresses  is  very 
important.  For  the  same  direction  of  shear  stress  the  two  notches  in  the  upstep  position 
give  a  result  close  to  the  one  obtained  with  one  notch;  indicating  very  little  coalescence  of 
the  damage  zones  of  the  two  notches.  On  the  other  hand  such  damage  zones  interfere  and 
coalesce  when  the  cracks  are  in  the  downstep  position,  leading  to  a  weakening  of  the 
specimen.  In  such  a  case,  for  a  given  distance  d,  e  plays  an  important  role  in  the  weakening 
process.  It  is  obvious  that,  as  d  increases,  the  behavior  tends  towards  that  of  the  single 
crack  specimen.  For  all  specimens  however,  it  seems  that  the  cracks  and  their  position  do 
not  affect  the  behavior  until  fifty  percent  of  the  stresses  have  been  applied.  Shear  bands 
become  visible  to  the  naked  eye  shortly  after.  The  corresponding  strain  is  about  0.5 
percent.  Photograph  taken  during  the  tests  and  showing  the  shear  bands  will  be  discussed 
in  Chapter  IV. 

3.7  Resonant  Column  Tests. 

A  number  of  resonant  column  tests  of  the  fixed-free  type  were  conducted  on  EPK 
clay  specimens  with  and  without  notches.  Some  were  conducted  on  specimens  with  two 
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notches.  Complete  description  and  results  of  those  tests  can  be  found  in  the  appendix. 
Young’s  modulii,  shear  modulii  and  damping  ratio  are  listed  for  each  test.  They  indicate 
that  for  the  levels  of  strain  involved  ( <  lO*^)  the  cracks  have  no  effect  on  the  response. 
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Normalized  Octahedral  Shear  Stress  vs.  Octahedral  Shear 
Strain;  EPK;  fi  =  0 
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Figure  3.1d  Normalized  Octahedral  Shear  Stress  vs.  Octahedral  Shear 
Strain;  EPK;  p  =  60° 
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Figure  3.2a  Normalized  Pore  Pressure  Developed  vs.  Octahedral  Shear 
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Figure  3.2c  Normalized  Pore  Pressure  Developed  vs.  Octahedral  Shear 
Strain;  EPK;  /3  =  45" 
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igure33a  Normalized  Octahedral  Shear  Stress  vs.  Octahedral  Shear 
Strain;  H121 
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Influence  of  the  Inclination  of  the  Major  Principal  S 
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igure3.5a  Comparison  of  Trimmed  and  Uniform  Specimens 
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Figure  3.6a  Bifurcation  and  Resulting  Shear  Bands;  0  =  30° 


Figure  3.6az  Bifurcation  and  Resulting  S 
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Figure  3.6b  Bifurcation  and  Resulting  Shear  Bands; 
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Comparison  Between  Notched  and  Unnotched  Specimens 
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Figure  3.8a  Comparison  Among  Various  Modes  of  Testing 
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igure  3.10a  Normalized  Octahedral  Shear  Stress  vs.  Octahedral  Shear 
Strain;  Overconsolidated  EPK 


Figure  3.10b  Normalized  Pore  Pressure  Developed  vs.  Octahedral  Shear 
Strain;  Overconsolidated  EPK 
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Figure  3.11a  Bifurcation  and  Resulting  Shear  Bands; 
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Figure  S.llbzBifurcation  and  Resulting  Shear  Bands;  /S  =  45* 
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Figure  3.11c  Bifurcation  and  Resulting  Shear  Bands;  p  - 


Figure  3.1  Icz  Bifurcation  and  Resulting  Shear  Bands; 
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Comparison  Between  Notched  and  Unnotched  Specimens 
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Figure  3.16  Decay  of  Shear  Modulus  and  Increase  in  Damping  of  an 
Anisotropic  EPK  with  an  Inclined  Crack 
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e  3.17a  Influence  of  the  Number  of  Cracks  and  Their  Relative  Position 
on  the  Strength 


e  3.17b  Influence  of  the  Number  of  Cracks  and  Their  Relative  Position 
on  the  Pore  Water  Pressure 


CHAPTER  IV 


KINEMATICS  OF  DEFORMATION,  CRACK  AND  DAMAGE  PROPAGATION 

Photographs  taken  during  and  at  the  end  of  many  of  the  tests  listed  in  Tables  2.1  to 
2.9  allow  one  to  study  the  pattern  of  crack  propagation,  shear  bands  development  and  the 
relative  displacement  of  the  various  parts  of  the  specimens.  One  of  the  difficulties  in 
developing  constitutive  equations  is  in  the  measurements  of  strains  within  the  soil 
specimens.  Once  a  shear  line  of  slip  band  is  formed,  whether  initiated  by  the  existence  of  a 
crack  or  not,  different  parts  of  the  specimen  under  test  move  as  blocks.  Bifurcation 
theories  and  models  have  been  developed  to  mathematically  describe,  indeed  predict  those 
phenomena  in  various  materials.  The  observations  described  in  this  chapter  are  but  a  first 
step  in  developing  a  model  suitable  to  clay  soils. 

Discontinuities  that  occur  during  deformation  have  been  classified  by  Morgenstem 
and  Tchalenko  (1967),  as  strain  discontinuities  and  displacement  discontinuities  (Fig.  4.1). 
Strain  discontinuities  generally  precede  the  displacement  ones  in  front  of  the  crack  tip  as 
well  as  in  front  of  the  propagating  shear  bands.  The  extensive  photographic  record  made 
during  deformation  and  failure  allows  us  to  follow  the  kinematics  of  the  specimen  while 
relating  it  to  the  resistance  of  the  material.  As  previously  stated,  this  record  was  digitized 
and  stored  in  our  data  base  so  that  a  pattern  of  slip  lines  can  be  directly  recalled  and 
plotted  without  referring  to  the  original  photograph. 

4.1  Specimens  without  Cracks  under  Static  Loading. 

Whether  it  is  under  axial  loads  or  under  a  combination  of  axial  and  torsional  loads, 
the  maximum  stress  is  reached  before  slip  lines  form.  Close  to  the  peak  the  deformations 
usually  occur  along  a  large  number  of  slip  surfaces  often  uniformly  distributed  along  the 
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length  of  the  cylinder.  Such  slip  surfaces  could  not  be  clearly  seen  in  the  past  because  of 
the  presence  of  filter  paper,  which  inhibited  the  formation  of  some  while  emphasizing 
others;  specially  where  the  paper  is  slit  to  decrease  its  resistance.  The  photogr^hs  in  the 
following  figures  illustrate  those  points: 

Fig.  4.2  shows  the  slip  lines  formed  on  a  specimen  of  H121  clay  subjected  to  torsion; 

Fig.  43a  shows  the  slip  lines  formed  on  a  specimen  of  EPK  clay  with  which  the 
external  filter  was  made  of  vertical  paper  strips.  Notice  the  shear  bands  between  the  strips; 

Fig.  4.3b  shows  the  vertical  shear  bands  after  very  large  torsional  deformations  in  an 
EPK  clay  specimen;  the  localization  being  along  the  once  vertical  filter  strips. 

Fig.  4.4  shows  the  slip  bands  that  develop  where  the  filter  paper  is  slit  horizontally 
to  decrease  its  resistance  to  torsion. 

Figs.  3.6  which  were  partially  commented  on  in  Chapter  HI  show  the  slip  bands  for  a 
normally  consolidated  clay  under  two  different  stress  paths  and  Figs.  3.11  show  the  slip 
bands  for  an  overconsolidated  clay  under  three  different  stress  paths.  No  outer  filter  paper 
was  used  and  the  grid  brings  out  the  discontinuities  that  take  place.  All  stress  paths  used  in 
those  tests  resulted  in  the  kind  of  discontinuities  illustrated  in  Fig.  4.1.  The  photographs  in 
those  two  sets  of  figures  were  taken  at  the  end  of  the  test;  in  other  words  after  the  release 
of  the  stresses  acting  on  the  specimens.  Some  rebound  obviously  occurred  but  one  would 
expect  that  all  parts  would  have  suffered  some  permanent  displacement.  Figs.  3.6  which 
correspond  to  a  normally  consolidated  clay  support  this  expectation.  However  Figs.  3.11 
which  correspond  to  an  overconsolidated  clay  show  regions  between  the  slip  lines  where 
the  vertical  grid  has  very  little  if  no  inclination.  Also,  it  has  been  noticed  during  the  tests 
that,  at  the  instant  of  formation  of  the  shear  bands,  some  rebound  takes  place  in  the  rest  of 
the  specimen. 
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The  inclination  of  the  observed  shear  bands  when  compared  to  the  inclination  of  the 
principal  stresses  deserves  comments.  Table  4.1  lists  normally  consolidated  clays  from 
Tables  2.1  and  2.2  and  overconsolidated  clays  from  Table  2.4  together  with  their  value  of  p, 
the  inclination  of  the  shear  bands  with  the  horizontal  S,  and  e  the  inclination  of  the  shear 
bands  on  the  major  principal  stress.  6  is  seen  to  be  nearly  constant.  Fig.  4.5  shows  6 
plotted  versus  b  (b  =  sinj9).  Similar  results  were  obtained  by  Arthur  and  Dunstan  (1982) 
for  non-cohesive  soils.  This  observation  shows  that  the  dilatation  angle  can  be  neglected; 
which  is  to  be  expected  globally  when  the  test  is  undrained. 

The  formation  of  the  shear  bands  close  to  the  peak  stresses,  and  the  subsequent 
occurrence  of  a  deformation  made  of  relative  motion  of  blocks  along  those  bands,  forces 
one  question  the  validity  of  the  displacements  measurements  and  the  subsequent 
calculations  of  the  strains.  The  jumps  that  one  notices  in  the  grids  are  clear  indications  of 
the  total  lack  of  uniformity  at  and  beyond  the  peaks.  Measurements  during  the  strain 
softening  period  caimot  be  used  to  support  or  discredit  many  of  the  constitutive  equations 
attempting  to  model  the  behavior  of  the  material  dining  this  stage. 

Additional  photographs  of  shear  bands  and  their  digitized  plots  can  be  found  in  the 
following  sections. 

4  J  Specimens  with  One  Crack  under  Static  Loadings. 

For  an  angle  fi  =  0,  i.e.  for  specimens  in  triaxial  compression,  the  crack  hardly 
affects  the  stress-strain  behavior.  Combination  of  axial  and  torsional  stresses  however, 
result  in  a  shear  band  starting  at  and  propagating  from  the  tip  of  the  crack,  before  the  peak 
is  reached.  While  in  the  absence  of  cracks  we  have  a  very  large  number  of  shear  lines  all 
around  the  specimen,  here  the  number  is  limited  to  two  or  three  all  starting  at  the  crack  or 
emanating  from  each  other.  The  first  band  starts  at  the  crack  tip  with  a  slope  that  is 
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dictated  by  the  stress  concentration  and  the  directions  of  the  principal  stresses  there.  It  is 
the  steepest,  while  the  second  and  third  start  at  a  smaller  slope.  As  they  propagate 
however,  they  tend  to  run  parallel  to  one  another  as  the  influence  of  the  crack  vanishes 
with  distance. 

Figures  4.6a,  b  and  c  show  the  shear  bands  for  the  three  normally  consolidated 
tests  in  Table  2.3  corresponding  to  /3  =  30, 45  and  60  degrees  respectively.  Notice  the  very 
large  displacements  in  4.6a  and  4.6c.  This  is  due  to  the  fact  that  this  was  a  stress  controlled 
test  where  failure  occurs  in  a  catastrophic  way  if  left  unchecked.  Compare  now  the 
patterns  and  directions  of  the  slip  lines  obtained  in  Figs.  4.6  to  those  obtained  if  Figs.  3.6. 
While  there  seems  to  be  an  infinite  (very  large)  number  of  shear  bands  in  Figs.  3.6,  the 
presence  of  the  crack  initiates  and  results  in  one  or  two  main  shear  bands  which  propagate 
all  around  the  specimen  to  cause  failure.  The  fact  that  blocks  of  specimens  move  with 
respect  to  each  others  once  the  shear  bands  form  is  very  striking  when  the  crack  is  the 
initiator.  In  Figs.  4.6  many  of  the  vertical  lines  remained  nearly  vertical  after  the  load  was 
removed  and  do  not  show  discontinuities.  Also,  once  one  moves  away  from  the  crack  the 
inclination  of  the  shear  bands  is  nearly  the  same  whether  the  specimen  has  a  crack  or  not. 
This  is  a  clear  indication  that  while  the  stress  concentrations  near  the  tip  of  the  crack 
control  the  moment  of  initiation  and  the  initial  direction  of  the  shear  band,  their  influence 
quickly  vanishes  leaving  the  plasticity  of  the  material  in  command  of  the  kinematics  of  the 
specimen. 

In  Figs.  4.6  the  main  system  of  shear  bands  appears  to  be  made  of  a  primary  band 
which  develops  first  followed  by  secondary  bands.  In  the  figures,  we  shall  use  the  notation 
of  Morgenstem  and  Tchalenko  (1967),  in  referring  to  those  bands.  The  first,  or  primary, 
will  be  called  Sj,  followed  by  S2,  S3 ...  etc.,  in  the  order  they  were  observed  to  form.  Those 
plots  follow  each  photograph  showing  shear  bands. 
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Table  4.2  lists  the  inclination  of  the  shear  bands  emanating  from  the  cracks.  This 
table  is  to  be  read  in  conjunction  with  the  digitized  plots  so  as  to  see  the  level 
approximation  involved  in  such  numbers.  The  table  also  lists  the  "Crack  Sliding 
Displacement,"  (C.S.D.)  for  the  various  tests. 

In  Fig.  4.6a,  and  others  corresponding  to  combinations  of  compression  and  torsion, 
one  notices  vertical  shear  bands.  It  is  believed  that  such  shear  bands  materialize  at  large 
deformations  as  a  release  mechanism  for  excessive  shearing  stresses. 

The  photographs  show  the  magnitude  of  the  slippage  along  each  of  the  main  shear 
bands.  The  slippage  varies  from  band  to  band  and  along  each  band  as  indicated  by  the 
grid.  It  was  often  observed  that  when  a  primary  shear  band  could  not  propagate  any 
further  because  of  boundary  restrictions,  secondary  bands  would  emanate  from  it.  This  is 
seen  in  Fig.  4.6a  and  its  digitized  plot. 

The  case  of  Fig.  4.6c  is  special  in  that  the  shear  bands  form  and  propagate 
essentially  along  the  original  direction  of  the  crack:  Primary  and  Secondary  bands  overlap 
and  meet. 

Figs.  4.7a  and  b  show  the  results  of  the  same  test  conducted  under  controlled  stress 
and  under  controlled  deformation  respectively,  but  with  fixed  length.  The  pattern  of 
deformation  is  practically  identical  including  the  verticality  of  the  grid  lines  below  the 
crack.  One  major  slip  line  develops  and  stops  after  a  while;  a  second  and  a  third  one  take 
over  until  failure  takes  place.  Notice  that  the  bands  are  more  inclined  than  those 
appearing  in  Fig.  4.6b  because  of  the  axial  stresses  that  develop  when  the  specimen  is 
prevented  from  increasing  in  length.  Once  the  main  shear  band  and  its  first  or  second 
offshoots  propagate,  practically  all  the  deformations  take  place  along  those  bands  with  the 
original  crack  as  part  of  the  slip  surfaces.  The  grid  shows  that  hardly  any  permanent 
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deformation  takes  place  in  the  rest  of  the  specimen.  Also  notice  that  Fig.  4.7a  indicates  a 
crack  shear  displacement  twice  as  large  as  that  seen  on  Fig.  4.7b. 

Inclined  cracks  results  in  shear  bands  that  start  at  the  tip  too  and  again,  at  a  short 
distance  from  the  crack,  take  a  direction  dictated  by  the  plastic  behavior  of  the  clay.  Fig. 
4.8a  shows  the  pattern  caused  by  a  downstep  crack  inclined  at  45  degrees  when  the 
specimen  is  subjected  to  torsion  and  free  to  deform  axially.  Additional  shear  bands  having 
no  physical  connection  with  the  crack  develop  at  the  peak.  This  can  clearly  be  seen  on  the 
photograph  as  well  as  on  the  digitized  plot. 

Fig.  4.8b  shows  the  results  of  a  test  identical  to  the  one  above  but  with  the  crack  in 
an  upstep  position;  the  slip  lines  at  the  front  and  tail  of  the  cracks  have  the  same 
inclinations  in  both  cases.  In  the  first  case,  however,  the  crack  is  part  of  the  surface  along 
which  the  blocks  move,  while  in  the  second  case  it  is  not. 

Fig.  4.8c  shows  a  downstep  crack  inclined  at  22.5  degrees.  Here,  a  shear  band 
emanates  from  the  tip  and  keeps  nearly  the  same  direction  as  that  of  the  crack.  Also,  there 
is  appreciable  crack  sliding  displacement  (C.S.D.)  when  compared  to  the  two  previous 
cases.  Only  one  shear  band  passing  through  the  crack  was  noticed,  in  other  words  the 
specimen  became  two  blocks  moving  with  respect  to  each  other. 

The  pattern  of  shear  bands  in  Figs.  4.8  is  extremely  instructive  as  to  the  influence  of 
the  inclination  of  the  principal  stresses  on  the  direction  of  the  crack.  In  the  case  of  the 
downstep  crack  (Fig.  4.8a)  torsion  results  in  a  combination  of  tension  aaoss  the  face  of  the 
crack  and  compression  along  it.  This  results  in  a  weak  point  at  the  tip  of  the  crack  resulting 
in  an  initiation  of  a  shear  band  there.  On  the  other  hand,  in  the  case  of  the  upstep  crack 
(Fig.  4.8b)  the  directions  of  tension  and  compression  with  respect  to  the  crack  are  reversed 
resulting  in  compression  across  the  faces.  The  specimen  reacts  as  if  there  was  no  crack  and 
the  normal  development  of  shear  bands  takes  place;  one  of  which  is  seen  to  be  passing 
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through  the  tip  (Fig.  4.8b).  In  both  the  cases  above  the  magnitude  of  the  Crack  Sliding 
EMsplacement  is  negligible.  The  crack  remains  closed  due  to  hydrostatic  stresses  and  the 
shear  along  the  faces  are  negligible. 

For  all  the  specimens  with  a  crack  the  shear  bands  show  up  first  before  the  peak  is 
reached.  This,  generally  causes  a  more  or  less  pronounced  kink  in  the  stress  strain  curves. 

The  shear  bands  observed  for  overconsolidated  clays  are  in  general  more  sharply 
defined  than  the  ones  obtained  for  normally  consolidated  clays.  The  tests  listed  in 
Table  2.5  illustrate  this  fact  specially  the  ones  conducted  under  controlled  rotation 
conditions: 

Figs.  4.9a,  b  and  c  show  the  shear  bands  for  p  =  30, 45  and  60  degrees  respectively. 
The  specimens  are  free  to  deform  axially.  Note  in  Fig.  4.9a  the  limited  number  of  shear 
bands  starting  at  the  tip  of  the  crack  and  spreading  circumferentially  around  the  specimen. 
Fig.  3.11a  which  corresponds  to  the  same  conditions  without  crack  has  practically  the  whole 
specimen  covered  with  slip  lines.  Similar  situations  are  encountered  for  the  cases  of 
torsion  and  combinations  of  tensions  and  torsion  as  shown  in  Figs.  4.9b  and  c  on  one  hand, 
and  Figs.  3.6a,  3.6b  and  Figs.  3.11a  and  3.11b,  and  3.11c  on  the  other  hand. 

It  is  also  interesting  to  notice  that  isotropically  overconsolidated  specimens  result  in 
a  diffuse  zone  of  shear  bands  starting  at  the  tip  of  the  crack,  but  with  limited  propagation. 
There  were  no  displacement  discontinuities  similar  to  the  ones  occurring  in  the 
consolidated  cases.  Fig.  4.9d  shows  this  pattern.  It  appears  that  random  orientation  of 
clusters  and  particles  impedes  the  formation  of  continuous  shear  bands. 

Figs.  4.10a  to  f  show  the  shear  bands  when  the  length  of  the  specimen  is  kept 
constant  and  when  the  size  of  the  crack  varies  from  a/1  =  0  to  a/1  =  0.4  (Table  2.5).  Fig. 
4.10a,  for  a/1  =  0,  shows  only  one  shear  line.  On  the  other  hand.  Fig.  3.11b  which 
corresponds  to  the  same  test  but  with  free  axial  deformation  indicates  the  presence  of  a 
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large  number  of  slip  lines.  It  appears  that  the  compressive  forces  that  are  axially  generated 
to  prevent  the  specimen  from  elongating  are  sufGcient  to  keep  the  multitude  of  slip  lines 
seen  on  Fig.  3.11b  to  form  and  become  active.  The  deformation  then  amounts  to  two 
blocks  slipping  on  each  others.  The  vertical  lines  of  the  grid  remain  continuous  except 
where  they  cross  the  slip  line.  By  looking  at  all  of  Figures  4.10a  to  4.10f  one  notices  that  as 
one  moves  away  from  the  crack  the  primary  shear  bands  adopt  the  inclination  of  the  one 
seen  in  4.10a  indicating  once  more  a  shear  behavior  dictated  by  the  plasticity  of  the 
material  and  not  by  the  stress  concentration  in  the  close  vicinity  of  the  crack.  It  is  to  be 
noticed  however  that  beyond  a/1  =  0.2  the  pattern  of  shear  bands  is  very  much  influenced 
by  the  size  of  the  cracks  with  respect  to  the  overall  size  of  the  sample.  It  is  not  known  how 
the  stress  distribution  in  the  specimen  is  affected  by  the  fact  that  we  do  not  really  have  a 
hollow  solid  cylinder  subject  to  pure  torsion  and  that  averaging  the  stresses  around  the 
circumference  is  not  a  valid  approximation. 

4J  Specimens  with  One  Crack  under  Dynamic  Loading. 

Slow  cyclic  loading  in  torsion  results  in  a  horizontal  crack  propagating  horizontally, 
preceded  and  surrounded  by  a  damage  zone.  The  slow  cyclic  loading  allows  the 
development  of  a  herringbone  typ>e  of  shear  bands  that  propagate  on  both  sides  of  the 
crack  as  it  advances.  It  is  essentially  the  Hrst  band  that  is  formed  in  the  static  test  that 
keeps  being  generated  as  the  crack  advances  and  the  stresses  are  reversed.  It  is  interesting 
to  note  that  most  of  the  damage  and  deformation  are  localized  in  a  strip  on  both  sides  of 
the  propagating  crack.  While  it  was  thought  that  this  kind  of  damage  could  only  occur  in 
overconsolidated  clays,  tests  conducted  on  normally  consolidated  clays  show  that  the  same 
characteristics  are  present  there  too;  even  though  the  magnitude  of  the  torque  was  about 
60  percent  of  the  peak  reached  under  static  loading.  (Remember  that  in  the  static  case  the 
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shear  bands  appear  at  about  80  percent  of  the  peak).  Figs.  4.11a  and  b  show  the  results  of 
two  tests  listed  in  Table  2.6  and  2.7,  both  at  fixed  length,  the  first  on  a  normally 
consolidated  specimen  and  the  second  on  an  overconsolidated  one.  Fig.  4.12  shows  the 
results  obtained  with  a  normally  consolidated  H121  clay  also  subjected  to  cyclic  torsion 
with  fixed  length.  All  cases  show  a  similar  herringbone  pattern  and  a  crack  propagating  in 
the  horizontal  direction.  The  pattern  obtained  when  the  specimen  is  fi'ee  to  deform  axially 
while  being  subjected  to  cyclic  torsion  is  not  any  different  as  indicated  by  Fig.  4.13.  Tests 
conducted  on  isotropic  overconsolidated  clays  led  to  a  damage  strip  significantly  narrower 
than  the  ones  obtained  with  K^,  overconsolidated  anisotropic  clays.  It  appears,  as 
previously  noticed,  that  randomly  oriented  clusters  do  not  result  in  well  defined  shear 
bands.  Fig.  4.14  clearly  illustrates  those  conditions. 

Cyclic  loading  on  a  specimen  with  an  inclined  crack  resulted  in  very  small 
herringbone  shear  bands  starting  at  the  tip  of  the  crack  and  spreading  in  a  horizontal 
direction;  showing  again  that  the  shear  bands  and  damage  zones  are  primarily  controlled  by 
the  system  of  stresses  and  the  plasticity  of  the  clay  (Fig.  4.15). 

Sometimes,  in  the  case  of  a  very  brittle  highly  overconsolidated  clay  under  small 
confining  pressure,  the  specimen  fails  under  a  single  cycle  showing  a  pattern  very  much  in 
line  with  the  predictions  of  classical  fracture  mechanics.  Fig.  4.16  is  a  photograph  of  a 
specimen  that  failed  under  those  conditions.  Notice  the  symmetry  of  the  fracture  lines  due 
to  tensile  stresses  that  develop  at  the  tip. 

The  impulse  loading  which  is  a  fast  monotonic  loading  does  not  allow  for  any  shear 
band  to  seriously  develop.  As  shown  in  Fig.  4.17a,  which  was  taken  at  the  end  of  the  test, 
most  of  the  vertical  lines  remain  continuous  and  at  best  the  type  of  discontinuity  shown  in 
Fig.  4.1a  is  present  near  the  crack.  Away  from  the  crack  the  shearing  strains  are  well 
distributed  showing  no  influence  of  crack  or  block  relative  motion.  Fig.  4.17b  shows  a 
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series  of  photographs  taken  with  an  SP2000  high  speed  camera  at  03m  sec  intervals  during 
a  similar  impulse  test  (Lesser  1989).  The  material  and  configuration  are  the  same  as  the 
ones  used  in  Fig.  4.17a,  but  the  clay  here  is  overconsolidated. 

4.4  Specimens  with  TVo  Cracks. 

Here  the  relative  position  of  the  cracks  plays  an  important  role  in  the  development 
of  the  shear  bands.  The  interaction  of  the  cracks  depends,  not  only  on  their  relative 
position  but  also  on  the  direction  of  the  applied  torsional  shear  stress.  In  this  investigation 
the  shearing  stresses  were  always  applied  in  one  direction  and  the  position  of  the  cracks 
changed  as  shown  in  Table  2.8.  The  cracks  when  placed  in  a  downstep  position  result  in  a 
tensile  instability  (local  necking)  which  is  quite  localized  and  allows  the  specimen  to  still 
resist  additional  loading  without  failing.  This  result  of  the  superposition  of  two  zones  with 
high  tensile  stress  was  illustrated  by  Vallejo  (1989).  Actual  calculation  of  stress 
distributions  will  be  given  in  Chapter  V.  Fig.  4.18  shows  the  interference  of  two  zones  of 
high  tensile  stress  at  the  tips  of  two  cracks.  This  is  the  first  test  listed  in  Table  2.8  and  for 
which  e  =  0.0  and  d  =  0.6  cm.  The  interference  appears  as  a  dark  hole  in  the  figures. 
Indeed,  at  this  location  the  specimen  becomes  locally  much  thinner  to  the  point  that  the 
inner  and  outer  membranes  touch  each  other.  The  very  local  tensile  failure  takes  place 
even  though  the  material  is  under  general  hydrostatic  compression  while  being  subjected  to 
torsion.  Fig.  4.19,  on  the  other  hand  shows  two  systems  of  shear  bands  each  developing  and 
spreading  on  its  own  without  visible  interference  effects.  Here  the  cracks  are  in  the  upstep 
position. 

With  the  cracks  in  the  downstep  position  and  e  increasing  there  comes  a  distance  at 
which  the  superimposed  tensile  stresses  are  not  large  enough  to  overcome  the  compressive 
spherical  stresses.  Figs.  4.20, 4.21  and  4.22  show  the  pattern  of  shear  band  development  as 
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e  takes  successive  values  of  0.6  cm,  13  cm  and  234  cm.  At  e  =  234  cm  the  crack  tips  are 
far  enough  from  one  another  so  that  no  tensile  instability  takes  place.  In  all  downstep 
cases,  the  shear  bands  pattern  resembles  quite  closely  that  observed  for  large  single  cracks. 
Indeed,  once  the  necking  takes  place  the  two  cracks  behave  as  a  single  large  one. 

Table  4.3  lists  the  inclinations  of  the  shear  bands  and  the  C.S.D.’s  for  the  case  of  two 

cracks. 

While  examining  the  photographs  notice  that  the  crack  shear  displacement  is  often 
nearly  equal  to  the  unit  size  of  the  grid  (0.6  cm).  There  are  reference  grid  lines  on  each 
photograph  that  indicate  the  size  of  the  displacement. 

4.5  Shear  Bands  Propagation  and  Their  Measurements. 

The  series  of  photographs  taken  at  the  initiation  of  a  shear  band  (at  the  tip  of  a 
crack),  and  during  its  subsequent  propagation  allow  one  to  plot  relations  involving  the 
length  of  the  band,  the  angle  of  rotation  of  the  top  of  the  specimen  with  respect  to  its  fixed 
bottom,  the  amount  of  slip  measured  along  the  band,  etc....  All  the  tests  for  which  such 
graphs  were  generated  are  listed  in  Table.  2.5 

Fig.  4.23  shows  the  length  of  the  shear  band  versus  the  angle  of  rotation  for 
specimen  18-1-15-91.  Fig.  4.24  shows  the  amount  of  slip  measured  along  the  first  shear 
band  at  the  end  of  the  test.  This  is  the  kind  of  displacement  discontinuity  shown  in  Fig. 
4.1c.  Fig.  4.25  shows  the  amount  of  slip  measured  along  the  first  shear  band  for  various 
angles  of  rotation. 

Figs.  4.26  and  4.27  correspond  to  specimen  66-05-29-91,  tested  under  variable  length 
conditions.  These  two  figures  differ  from  their  counterparts  in  Figs.  4.23  and  4.24  in  that 
they  show  the  slip  along  the  second  shear  band. 
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Figs.  4.28  and  4.29  correspond  to  specimen  16-12*13-90,  tested  under  fixed  length 
conditions. 

Figs.  430  and  431  correspond  to  specimen  70-06-14-91,  tested  under  combinations 
of  tension  and  torsion. 

In  all  the  cases  observed,  the  formation  and  propagation  of  the  shear  bands  was 
stable.  The  length  of  the  shear  band  varies  linearly  with  the  angle  of  rotation  (Figs.  4.23, 
4.26, 4.28  and  4.30).  The  second  shear  band,  after  its  initiation  propagated  at  a  rate  faster 
than  that  of  the  first  one;  it  also  became  the  dominant  one. 

The  appearance  of  additional  shear  bands  beyond  the  first,  depends  on  the  amoimt 
of  rotation  induced  in  the  specimen,  the  end  effects,  the  presence  of  the  inclusion  formed 
by  the  teflon  strips,  and  the  fact  that  the  first  shear  band  may  not  be  propagating  in  the 
direction  of  the  original  crack  (see  Fig.  4.8c) 

The  slip  that  takes  place  along  each  shear  band  varies  with  the  distance  from  the  tip 
of  the  crack  (Figs.  4.24,  4.27,  4.29  and  4.31).  While  it  was  possible  to  approximate  this 
relation  with  a  straight  line  in  Figs.  4.24  and  4.25,  this  is  not  the  case  in  Figs.  4.27  and  4.28. 
The  appearance  of  secondary  branches  affects  the  magnitude  of  the  slip  and  leads  to  a 
more  complex  relation  between  the  two  variables.  At  the  point  of  branching,  it  was  noticed 
that  the  magnitude  of  the  slip  along  the  main  band  is  nearly  equal  to  the  sum  of  the 
magnitudes  along  each  of  the  branches. 

The  measurements  and  observations  above  are  the  ingredients  needed  to  test  the 
validity  of  the  various  models  that  have  been,  or  will  be  proposed  for  bifurcation  and  shear 
banding  resulting  from  localization.  The  level  of  accuracy  of  such  measurements,  although 
reasonably  limited  by  the  fact  that  observations  are  made  from  outside  the  cell,  is  more 
than  sufficient  to  evaluate  the  suitability  of  a  behavioral  model. 
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TABLE  4.1  -  INCUNAllON  OF  SHEAR  BANDS  FOR  VARIOUS 
INCLINATIONS  OF  MAJOR  PRINCIPAL  STRESSES;  NO  CRACK 
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TABLE  42  INCLINATION  OF  THE  SHEAR  BANDS  AND  CRACK  SLIDING 
DISPLACEMENT  (CS.D.);  ONE  CRACK 


REF. 

LEFT 

SIDE 

RIGHT  SIDE 

TABLE 

DESIGNATION 

C.S.D. 

LABEL 

ANGLE 

LABEL 

ANGLE 

(cm) 

o 

o 

51-04-5-91 

2.20 

Si 

32 

Si 

30 

S2 

35 

S2 

26 

6V-09-26-90 

0.64 

Si 

30 

Si 

28 

S2 

25 

S2 

22 

52-04-09-91 

121 

Si 

0 

Si 

6 

S2 

3 

S2 

5 

2.3 

3f-02-06-90 

1.91 

Si 

25 

Si 

35 

S2 

22 

S2 

32 

6f-06-21-90 

0.95 

Si 

28 

Si 

33 

• 

S2 

20 

S2 

28 

56-04-17-91 

0.16 

Si 

20 

S2 

17 

S2 

20 

57-04-19-91 

0.8 

Si 

19 

Si 

20 

S2 

17 

58-04-23-91 

0.89 

Si 

20 

Si 

20 

2.9 

6F-09-17-90 

0.64 

51 

52 

30 

24 

51 

52 

36 

32 

m 

71-09-06-90 

0.64 

Slu 

23 

Slu 

23 

Sii 

8 

Su 

9 

COMMENTS:  Specimen  7a-09-06-90  is  isotropic  and  results  in  a  diffuse  zone  of  shear 
bands.  Subscripts  u  and  1  refer  to  upper  and  lower  bounds  respectively. 


(continued ....) 
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TABLE  42 


REF. 

TABLE 

DESIGNATION 

C.S.D. 

(cm) 

LEFT 

SIDE 

RIGHT  SIDE 

LABEL 

LABEL 

ANGLE 

o 

68-06-05-91 

0.89 

Si 

38 

Si 

39 

S2 

33 

S2 

33 

66-05-29-91 

0.95 

Si 

24 

Si 

30 

S2 

25 

S2 

23 

70-06-14-91 

0.36 

Si 

6 

Si 

4 

S2 

-8 

S2 

-7 

16-12-13-90 

0.45 

Si 

30 

Si 

28 

S2 

21 

S2 

18 

2.5 

18-01-15-91 

0.43 

Si 

29 

Si 

35 

S2 

26 

S2 

33 

15-11-28-90 

0.64 

Si 

30 

Si 

34 

S2 

29 

S2 

34 

13-11-07-90 

0.71 

Si 

31 

Si 

37 

S2 

27 

S2 

31 

14-11-20-90 

0.76 

Si 

27 

Si 

33 

S2 

27 

S2 

27 

TABLE  43  INCLINATION  OF  THE  SHEAR  BANDS  AND  CRACK  SLIDING 
DISPLACEMENT  (CS.D.);  TWO  CRACKS 


REF. 

TABLE 

DESIGNATION 

C.S.D. 

(cm) 

LEFT 

SIDE 

RIGHT  SIDE 

LABEL 

ANGLE 

o 

LABEL 

ANGLE 

o 

81-07-03-91 

0.51 

27 

r 

24 

Si 

33  _J 

S2 

28 

82-07-09-91 

0.25-0.50 

Si 

27 

Si 

26 

S2 

14  ^ 

S2 

18 

30 

^Sj 

22 

S2 

20 
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Types  of  Discontinuity  in  Defonnationf2l 


(a) 

(b) 

(c) 

(d) 
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Strain  Discontinuity  SS 
Displacement  Discontinuity  DD 

Magnification  of  Central  Portion  of  DD  Showing  Two 


Strain  Discontinuities  f^f^ 
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Figure  4.1  Morgenstem  and  Tchalenko’s  Definitions 
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Figure  4  Jb  Influence  of  the  Filter  Paper  on  the  Formation  of  Shear  Bands 
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Figure  4.4  Influence  of  the  Filter  Paper  on  the  Formation  of  Shear  Bands 
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Figure  4.5  Inclination  of  the  Shear  Bands 
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Figure  4.6a  Shear  Bands  for  One  Crack  and  p  =  30° 


Figure  4.6az  Shear  Bands  for  One  Crack  and  -  30° 
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Figure  4.6b  Shear  Bands  for  One  Crack  and  /9  =  45' 


Figure  4.6bz  Shear  Bands  for  One  Crack  and  =  45° 
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Figure  4.6c  Shear  Bands  for  One  Crack  and  /3  = 


Figure  4.7az  Shear  Bands  for  One  Crack  Under  Torsion  with  Fixed  Length 


Figure  4.7b  Shear  Bands  for  One  Crack  Under  Torsion  with  Fixed  Length 
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Figure  4.7bz  Shear  Bands  for  One  Crack  Under  Torsion  with  Fixed  Length 
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Figure  4.8az  Shear  Rands  for  an  Inclined  Crack  in  the  Downstep  Position; 
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Figure  4.8b  Shear  Bands  for  an  Inclined  Crack  in  the  Upstep  Position; 
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Figure  4.8bz  Shear  Bands  for  an  Inclined  Crack  in  the  Upstep  Position; 
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Figure  4.8c  Shear  Bands  for  an  Inclined  Crack  in  the  Downstep  Position; 


Figure  4.8cz  Shear  Bands  for  an  Inclined  Crack  in  the  Downstep  Position; 
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Figure  4.9a  Shear  Bands  for  One  Crack  aud  /3  =  30° 
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Figure  4.9b  Shear  Bands  for  One  Crack  and  /?  =  45' 


Figure  4.9bz  Shear  Bands  for  One  Crack  and  p  =  45° 
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Figure  4.9c  Shear  Bands  for  One  Crack  and  /3  = 


moo 

^  «H  « 
•  •H 
Q  h 


••  ••  •• 

O  H  ^ 

•p  o  .d 

•  M 
n  -p-H 
d  o 
O  H 

o 


Oi 

o 

U  m 
4i  I 
O  P$ 

n  o 
H  o 


o 

m 

d 

o 

o 


o  K 

Pi 

o 

Ot 

o 


••  •• 

O  H 

§d 
•H 
2  M 
9 

P  P 
CR  d 
O  2 


4.44 


Figure  4.9d  Shear  Bands  for  One  Crack  Under  Torsion  with  Fixed  Length; 
Isotropic  Clay 
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Figure  4.10a  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands;  a/I  =  0.0 


Figure  4.10az Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands;  a/I  =  0.0 
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Figure  4.10b  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands: 
a/I  =  0.10 
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Figure  4.10c  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands; 
a/I  =  0.18 


Figure  4.10cz  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands; 
a/1  =  0.18 
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Figure  4.10(1  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands: 
a/1  =  030 


Figure  4.10dzlnfluence  of  the  Size  of  the  Crack  on  the  Shear  Bands; 
a/I  =  OJO 


Figure  4.10e  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands; 
a/I  =  0 


Figure  4.10ez  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands; 
a/I  =  0J5 


e 


M  O 

H 


Q  h 


«  H  4i 
•M  043 
•  I4  » 
A 

0  • 

gH 


O  « 

S  o 


a 

o 

o 


O  K 
A  A 


M 


ii 


Z  M 

• 

4i  4J 
«  0 

g* 


Figure  4.10f  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands;  a/I  =  0.4 


Figure  4.10fz  Influence  of  the  Size  of  the  Crack  on  the  Shear  Bands;  a/I  =  0.4 
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Figure  4.11b  Crack  Propagation  and  Damage  Zone  for  a  K,,  Overconsolidated 
Clay  (EPK) 


Fiffure412  Crack  Propagation  and  Damage  Zone  for  a  Kq  Normally 
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Figure  4.13  Crack  Propagation  and  Damage  Zone  for  a  Overconsolidated 
Clay  (EPK) 
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Figure  4.14  Crack  Propagation  and  Damage  Zone  for  an  Isotropic 
Overconsolidated  Clay 
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Figure  4.15  Damage  Zone  for  an  Inclined  Crack 


Figure  4.16  Tensile  Failure  under  Cyclic  Torsional  Load 
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Figure  4.17a  Shear  Bands  and  Distortion  under  Impulse  Loading 
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Figure  4.18  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Downstep 
Position;  e  =  0 


Figure  4.18z  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Downstep 
Position;  e  =  0 
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Figure  4.19  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Upstep 
Position;  e  =  0 
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Figure  4.19z  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Upsten 
Position;  e  =  0 
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Figure  4.20  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Upstep 
Position;  e  =  0.6cm 
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Figure  4.21  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Upstep 
Position;  e  =  Uem 


Figure  4.21z  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Upstep 
Position;  e  =  Ucm 
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Figure  4.22  Shear  Bands  and  Interaction  of  Two  Cracks  in  the  Upstep 
Position;  e  =  2.54cni 
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Figure  4.22z  Shear  Bands  and  Interaction  of  TSvo  Cracks  in  the  Upstep 
Position;  e  =  2.54cni 
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Figure  4.23  Length  of  First  Shear  Band  for  Various  Levels  of  Rotation 
(18-1-15-91) 


Slip  measured  along  Shear  Band 
at  the  end  of  the  test 
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Figure  4.24  Amount  of  Slip  Along  the  First  Shear  Band  (18-1-15-91) 


Slip  along  the  first  Shear  Band 
for  various  levels  of  Rotation 
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Figure  4.25  Slip  Along  the  First  Shear  Band  for  Various  Levels  of  Rotation 
(18-1-15-91) 
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Figure  4.26  Length  of  Shear  Bands  for  Various  Levels  of  Rotation 
(66-5-29-91) 
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Figure  4.27  Amount  of  Slip  Along  the  Shear  Bands  (66-5-29-91) 
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Figure  4.28  Length  of  Shear  Bands  for  Various  Levels  of  Rotation 
(16-12-13-90) 
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Figure  4.29  Amount  of  Slip  Along 
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Figure  4 JO  Length  of  First  Shear  Band  for  Various  Levels  of  Rotation 
(70-6-14-91) 


Slip  measured  along  Shear  Band 
at  the  end  of  the  test 


Figure  4J1  Slip  Along  the  First  Shear  Band  for  Various  Levels  of  Rotation 
(70-6-14-91) 


CHAPTERS 


APPLICATION  OF  FRACTURE  MECHANICS  AND  PLASTICITY  THEORY 


The  observations  made  in  the  last  two  chapters  indicate  that  neither  fracture 
mechanics  theories  alone  nor  plasticity  theories  alone  can  describe  the  complete 
phenomenon.  This  was  pointed  out  by  Bazant  and  Mazars  (1990)  who  suggest  that  the 
problem  was  one  of  localization  of  damage  and  stability.  It  is  instructive  however  to  start 
with  the  application  of  linear  fracture  mechanics  and  proceed  from  there  to  combining 
fracture  mechanics  and  plasticity  to  the  material  at  hand. 

Using  the  results  obtained  from  the  tests  conducted  on  specimens  with  cracks  of 
various  a/1  ratios  (Table  5)  one  can  follow  the  standard  procedures  of  linear  fracture 
mechanics.  Fig.  5.1  shows  the  torque  versus  the  angle  of  rotation  for  a/1  varying  from  0  to 
0.4.  The  work  done  is  shown  in  Fig.  5.2  where  the  curves  are  seen  to  be  quite  close  to  each 
others  specially  at  small  angles  of  rotation.  For  ratios  a/1  >  0.1  all  the  curves  seem  to 
coincide  indicating  that  the  size  of  the  crack  has  little  effect  on  the  energy.  The  J  integral  is 
shown  in  Fig.  5.3  and  reflects  the  results  of  the  two  previous  figures. 

There  is  no  doubt  that  the  plasticity  of  the  material  plays  an  important  role.  Once 
the  shear  bands  start  forming  and  reach  some  length  they  control  the  behavior  with  little 
influence  from  the  original  size  of  the  crack.  However,  stress  distributions  calculated  using 
linear  elasticity  can  be  taken  advantage  of  to  explore  zones  of  high  stress  concentrations  in 
the  vicinity  of  the  crack’s  tip. 


5.1 


5.1  Stress  Concentration  Around  Cracks  Tips. 

Qassical  solutions  in  tenns  of  stress  intensity  factors  involve  ^proximations  that 
restrict  their  use  to  the  very  close  vicinity  of  the  crack.  However,  Muskhelishvili  (1953) 
gives  the  general  expressions  for  the  stress  distribution  in  an  infinite  plate  with  a  crack 
subjected  to  stresses  &  N2  at  infinity.  Using  a  standard  notation  of  elasticity  theory. 
(Fig.  5.4a)  we  have 
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<r  -  i<r  =  ^(2)  +  £1(5)  +  (2  -  2)^' (2) 
yy  xy 

where  the  complex  potential  ^(2)  and  £1(2)  are  given  by: 
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^(z)  *  - — IZZ - 
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2v4^^  2 


r  =  -  (N  +  N  ) 

4  '  2 


r'  =  -  -  (N  -  N  )  e 
2  '  " 


Thus, 
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<r  -  icr 
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^  -  a 
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For  the  shear  conditions  shown  in  Fig.  5.4b. 


r  =  0  F'  =»  iT.  f'  =  -iT  5.5 


Thus,  Eqs.  5.4  become 

iz 

<r  *  <r  =  -2t  Re  - 

XX  yy 


ia^z  +  iz  (2z^  -  3a^) 

0*  -  0-  +2i(r  *  T  -  5.6 

yy  XX  xy  ,2  2.3/2 

(z  -  a  ) 


For  a  combination  of  axial  lateral  and  shearing  stresses  superposition  is  used.  The 
crack  is  assumed  to  transfer  normal  stresses  and  stress  distributions  due  to  compressions 
are  not  affected  by  the  presence  of  the  crack.  Fig.  5.5  illustrate  the  superposition  used  in 
the  calculations.  The  third  dimension  is  ignored. 

An  appropriate  software  was  used  to  calculate  the  distribution  of  stresses  and  plot 
contour  lines  of  equal  shearing  stress  normal  stress  Oj^,  normal  stress  Oyy,  minimum 
principal  stress  03,  maximum  principal  stress  a^,  maximum  shear  stress  Otmax  equal 
mean  stress  Some  3D  graphs  were  also  made  to  show  how  fast  the  effect  of  the  stress 
concentration  at  the  tip  of  the  crack  died  down.  (This  effect  will  decrease  somewhat  for 
elastic-plastic  analysis.)  The  3D  graphs  also  allow  one  to  estimate  the  applicability  of 
equations  valid  for  the  infinite  medium  to  the  case  of  the  hollow  cylinder. 

Figs.  5.6a  to  5.6o  refer  to  a  pure  torsion  test  where  the  applied  shear  stress  is  equal 
to  1  kPa.  On  those  graphs  the  length  of  the  half  crack  is  2.54cm  (2  in.) 

Figs.  5.7a  to  5.7o  refer  to  combinations  of  tensile  axial,  hydrostatic  and  torsional 
stresses,  namely  to  Test  No.  52-4-9-91  in  Table  3.  The  applied  torsional  shear  stress  is 
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217.1430  kPa  and  the  applied  normal  stresses  are  192.441  kPa  axially  and  375.67  kPa 
laterally. 

Figs.  S.Sa  to  S.8o  refer  to  a  pure  torsion  test  under  hydrostatic  conditions,  namely  to 
test  No.  6V-9-26-90  in  Table  3.  The  applied  torsional  shear  stress  is  243225  kPa  and  the 
applied  normal  stresses  are  336.03  kPa.  This  last  pressure  of  course  is  the  effective 
spherical  stress  acting  on  the  specimen. 

Figs.  5.9a  to  5.9o  refer  to  combinations  of  compression  axial,  hydrostatic  and 
torsional  stresses,  namely  to  test  No.  51-4-5-91  in  Table  3.  The  applied  torsional  shear 
stress  is  275.17  kPa  and  the  applied  normal  stresses  are  543.86  kPa  axially  and  302.9  kPa 
laterally. 

Many  of  the  graphs  in  Figs.  5.7  to  5.9  indicate  zones  of  tension  in  spite  of  the 
reasonably  high  hydrostatic  stress  acting  on  the  specimen.  Such  zones,  as  well  as  the  ones 
exhibiting  high  shearing  stresses  would  be  prime  areas  in  which  the  crack  or  the  shearing 
bands  would  propagate.  The  graphs  showing  the  contour  lines  for  maximum  shearing 
stresses  should  give  clear  indications  on  the  direction  of  r*  shear  band,  should  one  appear  at 
the  tip  of  the  crack;  if  the  material  were  a  Tresca  or  a  Mises  material  in  the  plastic  range. 
Indeed  such  directions  are  shown  on  Figs.  5.7k,  5.8k  and  5.9k  by  a  slashed  line.  The 
directions  of  the  shear  bands  observed  in  the  last  chapter  and  shown  in  Figs.  4.6c,  b  and  a 
show  that  this  is  not  the  case  and  point  to  the  need  for  an  elasto-plastic  analysis  taking  into 
account  the  frictional  nature  of  the  material. 
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5J2  Finite  Elements  Elasto>Plastic  Analysis. 

The  finite  element  code  ABAQUS  allows  one  to  study  the  stress  distribution  and 
the  propagation  of  the  plastic  zone  as  stresses  on  the  specimen  increase.  Fig.  5.10  shows 
the  size  of  the  plate  (developed  hollow  cylinder)  for  which  the  calculations  were  made. 
Notice  that  the  crack  has  a  blunted  end  to  represent  the  thickness  of  the  Teflon  inserts. 
The  proper  size  of  mesh  was  obtained  by  calculating  J  integrals  for  three  different  mesh 
configurations  and  judging  the  change  when  passing  from  one  mesh  to  the  other.  A  mesh 
with  eight  nodes  isoparametric  elements  was  chosen.  A  total  of  480  elements  with  1626 
nodes  was  used. 

For  the  elasto-plastic  analysis  a  Young’s  modulus  of  20.7  x  10^  kPa  and  a  Poisson 
ratio  of  0.3  were  used  together  with  either  a  Mises  or  a  Drucker-Prager  (DP)  yield 
criterion.  The  angle  for  the  Drucker-Prager  criterion  is  given  by  (Fig.  5.11) 

6  sin 

tanp  =  - 

3  -  sin 

where  0 '  is  the  effective  angle  of  friction  chosen  as  33.3  degrees.  The  value  of 

6  cos 

d  =  c  - 

3  -  sin 

is  chosen  to  correspond  to  a  cohesion  of  172.5  kPa. 

Figs.  5.12  and  5.13  show  the  propagation  of  the  yield  zone  for  the  Mises  and  DP 
criteria  respectively.  The  test  is  a  pure  torsion  test.  Notice  the  direction  of  propagation 
when  one  changes  from  a  non-frictional  to  a  frictional  material.  As  expected  the  Mises 
criterion  gives  a  direction  of  propagation  that  is  along  that  of  the  crack.  On  the  other  hand 
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the  frictional  material  shows  a  definite  inclination  of  the  direction  of  the  crack.  This  is 
what  was  observed  in  our  experiments.  Therefore,  the  DP  criterion  will  be  used  in 
analyzing  results  obtained  from  tests  imder  combined  stresses. 

The  first  step  is  to  compare  the  elastic  solution  with  the  elasto-plastic  one.  The 
results  of  the  elastic  solution  for  the  pure  toision  cases  are  shown  in  Figs.  S.14a,  b  and  c. 
[They  compare  exceedingly  well  with  the  Muskhelishvili  solutions  illustrated  in  Fig.  5.6g,  i 
and  k].  The  results  of  the  elasto-plastic  solution  are  shown  in  Figs.  5.15a,  b  and  c.  Both  sets 
of  figures  give  lines  of  equal  aj,  03  and  oj  -  03,  respectively.  On  notices,  particularly  on 
the  graphs  representing  that  the  tension  zone  is  much  smaller  for  the  elasto-plastic 
analysis  than  it  is  for  the  elastic  analysis  (see  contour  No.  29  along  which  there  is  tension). 
At  a  short  distance  from  the  tip,  however,  the  two  sets  of  contour  lines  are  quite  similar  for 
all  cases. 

A  better  appreciation  of  the  differences  between  the  two  analyses  can  be  obtained 
from  looking  at  Figures  5.16a  and  b  and  Figures  5.17  a  and  b.  Figure  5.16a  and  b  show  the 
values  of  oi  -  03,  oj  and  03  along  a  horizontal  line  starting  at  the  tip  of  the  crack  for  the 
elastic  and  the  elasto-plastic  analyses  respectively.  Figs.  5.17a  and  b  show  the  same 
quantities  along  a  vertical  line  at  the  tip  of  the  crack.  One  notices  that  the  distribution  in 
the  vertical  direction  is  more  affected  by  the  type  of  analysis  than  it  is  along  the  horizontal 
direction.  Also  seems  to  be  the  quantity  that  is  most  affected  in  both  cases.  It  is 
compressive  near  the  tip  of  the  crack  in  Fig.  5.17b  and  tensile  in  5.17a;  in  spite  of  the  high 
hydrostatic  pressure  of  nearly  340  kPa  acting  on  the  material.  The  plastic  zone,  which  very 
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quickly  appears  at  the  tip  of  the  crack,  eliminates  the  high  stress  concentration  that  is  seen 
in  Fig.  5.17a. 

Three  tests  involving  combined  stresses  were  analyzed  elasto-plasticly  using  the  DP 
criterion.  A  compression-torsion  test  (51-4-5-91),  a  torsion  test  (6v-9-26-90)  and  a  tension 
torsion  test  (52-4-9-91);  aU  listed  in  Table  3. 

Figs.  5.18a,  b  and  c  show  the  spread  of  the  plastic  zone  for  the  compression-torsion 
test,  the  pure  torsion  test  and  the  tension-torsion  test.  The  loading  was  radial  and  was 
increased  to  a  value  of  shearing  stress  of  about  170  kPa  with  the  axial  load  following 
proportionally.  The  contour  lines  are  lines  of  equal  equivalent  plastic  strain  defined  as 

p  r/^de  dc 
e  =  V  3  ij  ij 
eq  J 

The  evolution  of  the  plastic  zones  suggest  that  when  the  material  localizes  it  will  result  in 
shear  bands  having  directions  similar  to  the  ones  observed  in  Chapter  4  i.e.  approximately 
35°,  22°  and  0°  for  three  tests.  The  high  concentration  of  contours  in  the  upper  part  of  the 
crack  tip  is  due  to  the  high  tension  existing  in  this  region. 

For  the  sake  of  completeness  the  contours  giving  aj,  03  and  oj  -  03  are  also  given 
in  Fig.  5.19,  520  and  5.21  for  the  three  tests.  In  those  figures  a,  b  and  c  refer  to 
compression-torsion,  pure  torsion  and  tension-torsion  respectively. 

Tables  5.1  and  52  allow  one  to  obtain  the  magnitudes  of  the  stresses  on  the  contour 
lines  shown  on  Figs.  5.12  to  521. 
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TABLE  5.1  STRESS  VALUES  ON  VARIOUS  CONTOUR  LINES 


no.  No. 

STRESS  ON 
CONTOUR  No.l 
kPa 

STRESS  ON 
CONTOUR  No.  30 
kPa 

INCREMENT 

kPa 

5.14a 

5.15a 

5.19a,  b,  c 

1378.0 

68.9 

45.14 

5.14b 

5.15b 

5.20a,  b,  c 

-1378.0 

68.9 

49.89 

5.14c 

5.15c 

5.21a,  b,  c 

-689.0 

344.5 

35.64 

Remark:  The  stress  on  any  contour  line  between  contours  No.  2  and  No.  30  is  calculated 
as  follows: 

Stress  =  (Contour  Number-1)  x  Increment  +  Stress  on  Contour  No.  1 
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TABLE  5.2  EQUIVALENT  PLASTIC  STRAIN  VALUES  ON 
VARIOUS  CONTOUR  LINES 


no.  No. 

STRAIN  ON 
CONTOUR  No.  1 

STRAIN  ON 
CONTOUR  No.2 

STRAIN  ON 
CONTOUR  No.  30 

INCREMENT 

5.12 

5.13 

5.18a,  b,  c 

3.44x10-^ 

3.44x10^ 

10^ 

3.45  X  10-^ 

Remark:  The  equivalent  plastic  strain  on  any  contour  line  between  contours  No.  2  and  30 
is  calculated  as  follows: 


Strain  =  (Contour  Number-2)  x  Increment  +  Strain  on  Contour  No.  2 
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Versus  Rotation 
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Figure  5.1  Influence  of  the  Size  of  the  Crack  on  the  Strength 


J-Integral  vs.  Rotation 
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Figure  53  J-Integral  vs.  Rotation 


Figure  5.4a  Notation 


Superposition  Model  for  Elastic  Analysis 
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Figure  5.5  Superposition  Model  for  Elastic  Analysis 


Contour  of  the  Shear  Stress  Oxy 
- Reference:  Unit  Shear  (pure— torsion) 


Shear  Stress 
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Figure  5.6b  3-D  Plot  of  Shear  Stress 


Contour  of  the  Normal  Stress  (jxx 
- Reference:  Unit  Shear  (pure— torsion) 
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Figure  5.6c  Lines  of  Equal  Normal  Stress 
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Figure  5.6d  3-D  Plot  of  Normal  Stress 
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Figure  5.6e  Lines  of  Equal  Normal  Stress 
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Figure  5.6f  S-D  Plot  of  Normal  Stress 


Contour  of  the  Maximum  Principal  Stress 

- Reference:  Unit  Shear  (pure-  torsion) 
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Figure  5.6g  Lines  of  Equal  Maximum  Principal  Stress 


Principal  Stress 


.6h  3-D  Plot  of  Maximum  Principal  Stress 
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Figure  5.6i  Lines  of  Equal  Minimum  Principal  Stress 


Figure  5.6j  3-D  Plot  of  Minimum  Principal  Stress 


Dlot  for  the  Maximum  Shear  Stress 
Reference:  Unit  Shear  (piire-torsion) 


e  5.61  3-D  Plot  of  Maximum  Shear  Stress 


Contour  of  the  Mean  Stress 
- Reference:  Unit  Shear  (pure—torsion) 
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Figure  5.6m  Lines  of  Equal  Mean  Stress 
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Figure 


Principal  Stresses  and  Their  Direction 

(Unit  Shear  Stress;  Pure— Torsion;  Scale:  1  unit  ==  10  kPa) 
-  Tension  -  Compression 
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Figure  5.6o  Principal  Stresses  and  Tlieir  Directions 


Stress  Calculation  for  EPK  #52:  4/9/91;  Tension-Torsion 
- Ko  N.C.“Hor.N0TCH*a/l»0.183“F.H.‘C.S.‘ecp  =  70-20  psi‘DTR60 

length  of  craclc  »  1.0000  inch 

initial_x-  -1.0000  stepx*  0.1000  finial_x»  2.0000 

initially*  1.0000  stepy*  -0.1000  finial_y»  -1.0000 

sxy  *  -217.1430  syy  =  -192.4410  Sxx  =  -374.6700  Seff  =  -374.6700 
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Contour  of  the 

-  EPK  #52-04-  ( 
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Figure  5.7a  Lines  of  Equal  Shea 
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Figure  5.7c  Lines  of  Equal  Normal  Stress 
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Figure  5.7d  3-D  Plot  of  Normal  Stress 
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Figure  5.7f  3-D  Plot  of  Normal  Stress 


Contour  of  the  Maximum  Principal  Stress 

- EPK  #52-04-09-91  (tension-torsion) 
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Figure  5.7g  Lines  of  Equal  Maximum  Principal  Stress 


Contour  of  the  Minimum  Principal  Stress 

- EPK  j5i52— 04— 09— 91  (tension— torsion) 
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Figure  5.7i  Lines  of  Equal  Minimum  Principal  Stress 
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Figure  5.71  3-0  Plot  of  Maximum  Shear  Stress 


Contour  of  the  Mean  Stress 

- EPK  #52-04-09-91  (tension-torsion) 


Mean  Stress 


plot  for  the  Mean  Stress 

EPK  #52-  04-09-91  (tension-torsion) 


5.44 


Figure  5.7n  3-D  Plot  of  Mean  Stress 


Principal  Stresses  and  Their  Direction 

(EPK  #52-04-09-91;  Tension-Torsion;  Scale:  1  unit  =  3500  kPa) 
-  Tension  -  Compression 
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Principal  Stresses  and  Their  Directions 


Stress  Calculation  for  EPK  #6v:  9/26/90;  Pure-Torsion 
- Ko  N.C.“Hor.N0TCH"a/l*0.183*F.H.*C.D.“ecp  *  70-20  psi"R45 

length  of  crack  »  1.0000  inch 

initial_x»  -1.0000  stepx*  0.1000  finial_x*  2.0000 

initial_y»  1.0000  stepy*  -0.1000  finial_y»  -1.0000 

sxy  *  -243.2250  Syy  =  -336.0300  Sxx  =  -336.0300  Seff  =  -336.0300 
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Contour  of  the  Shear  Stress 

- KPK  ^^6v-  09-  26-  90  (pure-  torsion) 
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Figure  5.8a  Lines  of  Equal  Shear  Stress 
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Figure  5.8b  3-D  Plot  of  Shear  Stress 


Contour  of  the  Normal  Stress 

- EPK  #6v— 09— 26— 90  (pure— torsion) 
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Figure  5.8c  Lines  of  Equal  Normal  Stress 
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Figure  5.8d  3>D  Plot  of  Normal  Stress 


Contour  of  the  Normal  Stress 

- EPK  #6v— 09— 26— 90  (pure— torsion) 
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Figure  5.8e  Lines  of  Equal  Normal  Stress 


Prindpa.  Stress 


Contour  of  the  Minimum  Principal  Stress 

- EPK  #6v— 09— 26—90  (pure— torsion) 


Minimum  Principal  Stress 


Figure  5.8j 


Figure  5.8k  Lines  of  Equal  Maximum  Shear  Stress 
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Figure  5.8m  Lines  of  Equal  Mean  Stress 


Principal  Stresses  and  Their  Direction 

(EPK  #6t-O9-26-0O:  Pure-Torslon;  Scale;  1  unit  =  4000  kPa) 
-  Tension  -  Compression 
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Figure  5.8o  Principal  Stresses  and  Their  Directions 


Stress  Calculation  for  EPK  #51;  4/5/91;  Compression-Torsion 
- Ko  N.C. “Hor.N0TCH“a/l=0. 183“F.H. *C.S. "ecp  =  70-20  psi*DCR30 

length  of  crack  =  1.0000  inch 

initial_x*  -1.0000  stepx=  0.1000  finial_x=  2.0000 

initial_y=  1.0000  stepy=  -0.1000  finial_y=  -1.0000 

Sxy  =  -275.1720  Syy  =  -543.8580  SXX  =  -302.9100  Seff  =  -302.9100 
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D  plot  for  Normal  Stress  0”xx 

EPK  #51—04—05—91  (compression— torsion) 
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Figure  5.9d  3-D  Plot  of  Normal  Stress 


oiitour  of  the  Maximum  Principal  Stress 

- EPK  ^^51— 04— 05— 91  (compression— torsion) 


Maximum  Principal  Stress 


plot  for  Minimum  principal  Stress 

EPK  #51—04—05—91  (compression— torsion) 
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Figure  5.9j  3-D  Plot  of  Minimum  Principal  Stress 


Coilloiir  ol'  til 

- El’K  //51-  04- 


5.73 


Figure  5.9k  Lines  of  Equal  Maxim 


3— D  plot  for  Mean  Stress 

EPK  #51—04—05—91  (compression— torsion) 
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Figure  5.9n  3-D  Plot  of  Mean  Stress 
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Figure  5.9o  Principal  Stresses  and  Their  Directions 
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Figure  5.11  Convention  for  the  Drucker-Prager  Criterion 
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Figure 


5.83 


Figure  5. 
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Figure  5. 


Figu 


Figure  5.15c  Lines  of  Equal  Principal  Stress  Difference;  Eiasto-Plastic 
Solution;  Pure-Torsion 
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e  5.16a  Elastic  Stress  Distribution  along  OA 
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Figure  5.16b  Elasto-Plastic  Stress  Distribution  along  OA 
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Figure  5.17a  Clastic  Stress  Distribution  along  O 
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Figure  5.20c  Lines  of  Equal  Principal  Stress  Difference;  DP  Model;  Pure- 
Torsion  Test 


Figure  5.21a  Lines  of  Equai  Mqjor  Principal  Stress;  DP  Model;  Tension- 
Torsion  Test 


Figure  5.21b  Lines  of  Equal  Minor  Principal  Stress;  DP  Model;  Tension- 
Torsion  Test 
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Figure  5.21c  Lines  of  Equal  Principal  Stress  Difference;  DP  Model;  Tension- 
Torsion  Test 


CHAPTER  VI 


FABRIC  CHANGE  AND  MICROMECHANISMS  OF  DEFORMATION 

6.1  Macroscopic  Observations. 

External  observation  of  shear  band  formation  during  static  tests  and  damage  zone 
propagation  during  dynamic  tests  resulted  in  a  clear  description  of  the  kinematics  of 
deformation  of  a  hollow  cylinder  subjected  to  various  systems  of  stresses.  For  internal 
examination  the  impregnation  technique  described  in  Chapter  n  was  used. 

Before  damage  can  be  characterized,  the  kaolinite  must  be  identified  in  its  natural 
state.  Clay  fabric  has  been  studied  among  others  by  Yong  and  Warkentin  (1975)  and 
Mitchell  (1976)  and  a  scale  hierarchy  in  the  microstructure  identified.  Initially  the  clay 
particles  (10*^m)  are  aggregated  or  flocculated  together  in  units  called  domains  (10*^).  The 
domains,  in  turn,  form  together  clusters  (lO'^m  to  lO’^m).  Similarly,  the  clusters  group 
together  to  form  peds  or  larger  clusters  (lO'^m).  The  large  clusters  are  large  enough  to  be 
seen  without  a  microscope. 

A  polished  section  from  an  untested  hollow  cylinder  impregnation  is  shown  in  Fig. 

6.1. 

In  this  figure  the  kaolinite  appears  white  and  the  impregnation  material  appear 
dark.  The  kaolinite  clusters  are  homogeneously  distributed  and  have  a  characteristic  size 
of  the  order  of  10"^  meters  which  compares  with  what  has  been  reported  by  others. 

Figs.  6.2  compares  two  radial  sections  cut  from  a  cylinder  of  EPK  clay  cyclicly 
loaded,  after  impregnation;  one  obtained  from  an  undamaged  region  and  the  other  in  firont 
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of  the  tip  of  the  crack.  The  undamaged  section  in  Fig.  6.2a  shows  the  previously  reported 
clusters  homogeneously  distributed  in  the  wax  matrix.  The  damaged  section  in  Fig.  6.2b 
indicates  substantial  morphological  changes  due  to  the  loading  condition.  In  this  region, 
one  notices  a  change  in  the  size  distribution  of  the  clusters.  A  few  coalesce  and  others  split 
into  smaller  units.  Also,  one  can  see  the  creation  of  horizontal  strips  along  which  the  clay 
clusters  have  agglomerated.  These  strip  constitute  surfaces  within  the  cylinder.  The 
process  producing  the  surfaces  has  been  called  by  Lesser  (1989)  "Localized  Strip 
Densification"  or  LSD  for  short.  Notice  the  interior  fabric  consisting  of  densified  strips 
intermixed  with  clusters. 

Similar  results  were  observed  with  the  specimens  of  H121  clay.  The  impregnation 
time  however  for  this  clay  was  quite  different  from  that  of  the  EPK. 

6,2  Strips  Characterization. 

The  first  step  in  the  characterization  process  involves  making  a  series  of  radial 
sections  through  a  damaged  zone  of  an  impregnated  cylinder.  The  sections  are  usually 
equally  spaced  at  approximately  5  mm  interval  and  start  near  the  notch  tip  and  continue 
through  the  entire  damage  zone.  Each  section  is  mounted  on  a  slide  and  polished.  It  is 
then  segmented  into  3mm  thick  bands  below  a  known  reference  plane  as  shown  in  Fig.  6.3. 
Each  slide  is  then  placed  under  the  stereoscope  and  the  image  analyzer  is  used  to  compute 
the  total  LSD  area  in  each  band.  These  areas  are  then  used  to  construct  a  histogram  for 
each  slide  as  illustrated  in  Fig.  6.4. 
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It  was  observed  that  the  LSD  area  distributions  are  generally  bell-shaped,  thus  a 
Gaussian  curve  is  used  to  approximate  the  distribution  in  the  axial  direction.  These 
histograms  and  their  approximating  distributions  are  illustrated  in  Fig.  6.5.  Lesser  (1989) 
gives  three-dimensional  plots  showing  the  evolution  of  the  LSD’s  as  a  function  of  the 
number  of  cycles. 

Figs.  6.6  show  the  clusters  and  LSD’s  in  an  isotropic  EPK  clay.  They  are  not  very 
different  from  the  ones  obtained  from  a  consolidated  clay. 

Figs.  6.7  show  undamaged  and  damaged  regions  of  anisotropic  H121  clay.  The 
undamaged  region  shows  clusters  flattened  by  one-dimensional  consolidation.  This  is 
different  from  what  was  noticed  in  the  case  of  EPK  clay.  In  the  damaged  zones,  the  LSD’s 
although  thicker  do  not  differ  significantly  from  their  counterparts  in  the  EPK  clay. 

Observations  in  the  surface  analyzer  made  on  radial  sections  of  impregnated 
specimens  show  no  formation  of  LSD’s  to  speak  of  in  static  tests.  The  high  degree  of 
reorientation  of  particles  and  clusters  in  cyclic  tests  is  primarily  due  to  the  significant 
amount  of  energy  to  be  dissipated  by  the  material. 

6  J  Microscopic  Observations. 

Impregnated  specimens  as  well  as  non-impregnated  ones  were  examined  in  an 
electron  microscope  in  the  laboratories  of  McGill  University  by  Dr.  R.  Yong  and  his  staff. 
It  was  found  that  the  impregnation  took  place  in  the  micropores.  In  addition,  horizontal 
and  vertical  sections  made  in  the  undamaged  parts  of  the  K^  consolidated  EPK  clay  and  of 
the  isotropic  EPK  clay  showed  no  preferential  arrangements  within  the  clusters.  The  clay 
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particles  were  extremely  small  and  randomly  oriented  there,  llie  arrangement  at  the 
cluster  level  appears  to  be  the  main  reason  for  the  anisotropic  mechanical  response  of  this 
material.  It  is  felt  that  the  model  illustrated  by  Yong  and  Warkentin  on  pages  86  to  88  of 
their  treatise  on  Soil  Properties  and  Behavior  applies  to  the  clay  at  hand  and  can  be  used  to 
explain,  at  least  partially,  the  behavior  in  the  damage  zone  in  front  of  the  crack  (Fig.  6.8). 

Figures  6.9  and  6.10  are  micrographs  showing  horizontal  and  vertical  sections  within 
the  clusters  in  the  over-consolidated  anisotropic  clay  prior  to  subjecting  it  to  cyclic  loading. 
No  preferential  directions  are  seen  indicating  an  isotropic  or  random  arrangement. 
Remember,  however,  that  when  tested  for  strength  this  clay  behaved  in  an  anisotropic  way. 

Hydrite  121  was  found  to  have  much  larger  particles.  Fig.  6.11  shows  that  this  clay 
has  tendencies  of  orientation  within  the  clusters  when  the  clay  is  one-dimensionally 
consolidated. 

6.4  Conclusions. 

The  observations  reported  in  this  chapter  indicate  that  the  deformations 
mechanisms  in  a  saturated  clay  occur  primarily  at  the  cluster  level.  The  appearance  of  the 
LSD’s  during  cyclic  loading  and  the  coalescence  of  the  cluster  show  the  potential  for 
rearrangement  at  the  level.  The  examinations  in  the  surface  analyzer  and  in  the  electron 
microscope  both  point  in  that  direction.  The  model  advanced  by  Yong  and  Warkentin 
seems  to  be  a  valid  one. 
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Figure  6.1  Impregnated  specimen  of  EPK 
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Figure  6^  Sections  of  Undamaged  and  Damaged  Kg  Overconsolidated  EPK 
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LSD  AREA 
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Figure  6.4  LSD  Histograms  at  Various  Locations  in  Front  of  the  Crack  Tip. 
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Figure  6.6  Sections  of  Undamaged  and  Damaged  Isotropic  EPK 
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Figure  6.7  Section  of  Undamaged  and  Damaged  Consolidated  11121 
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Figure  6.8  Clay  Fabric  Model 
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Figure  6.9  Photomicrograph  of  Horizontal  Section;  EPK 
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Figure  6.10  Photomicrograph  of  Vertical  Section;  EPK 


Figure  6.11  Photomicrograph  of  H121  (Yong  &  Warkentin) 
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CHAPTER  VII 


CONCLUSIONS 


This  research  has  endeavored  to  examine  the  behavior  of  saturated  clays,  in  a 
reasonably  complete  way,  bringing  together  elements  of  continuum  mechanics,  fracture 
mechanics  and  stability.  Localizations,  spontaneous  or  induced  by  the  presence  of  cracks, 
appear  as  shear  bands  along  which  most  of  the  deformations  take  place.  Under  cyclic 
loading  they  form  damage  zones  within  which  clusters  of  clay  suffer  serious  rearrangement. 

7.1  In  the  Area  of  Strength: 

It  was  found  that  for  consolidated  clays  the  response  to  stress  was  significantly 
influenced  by  the  inclination  of  the  principal  stresses  on  the  axis  of  symmetry  of  the 
material.  Bifurcations  in  the  form  of  shear  bands  invariably  appeared  near  the  peak  stress 
in  strength  experiments.  Initially  well-distributed  along  the  length  of  the  hollow  cylindrical 
specimens,  the  shear  bands  eventually  relinquish  their  share  of  the  deformation  to  one 
along  which  the  material  fails  in  a  catastrophic  way. 

The  presence  of  a  crack  (or  notch)  in  a  specimen  did  not  substantially  affect  its 
overall  strength.  Indeed,  stress-strain  curves,  for  specimens  with  or  without  cracks,  differed 
very  little  even  though  the  observed  pattern  of  displacement  was  far  fi'om  being  the  same. 
The  pore  water  pressures,  which  are  of  course  measured  at  the  specimens  ends  are  hardly 
affected.  One  notices  however  a  lack  of  smoothness  in  the  stress-strain  curves  related  to 
notched  specimens;  specially  when  the  crack  is  oriented  in  the  direction  of  the  forming 
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shear  bands.  Here  the  bands  start  forming  at  about  eighty  percent  of  the  peak  value  of  the 
stress. 

Under  cyclic  loading  anisotropic  clays  yield  an  unsymmetric  pattern  of  hysterisis 
loops,  whether  the  specimen  had  a  crack  or  not.  The  presence  of  an  inclined  crack 
however  caused  the  hysterisis  loop  to  adopt  a  symmetric  pattern  characteristic  of  isotropic 
clays. 

The  presence  of  two  cracks  in  close  proximity  to  each  other  led  to  interferences 
which  depended  on  their  relative  position  and  on  the  direction  of  the  applied  shearing 
stresses.  These  interferences  affected  the  value  of  the  peak  stress  as  well  as  the 
development  of  the  shear  bands.  The  bands  were  noticed  at  about  80  percent  of  the  peak 
stress. 

.  All  our  observations  led  us  to  conclude  that,  while  cracks  create  local  stress 
concentrations,  the  behavior  of  the  material  is  dictated  by  its  plasticity  which  becomes 
predominant  at  a  very  short  distance  from  the  tip. 

12  \n  the  Area  of  Bifurcation  and  Shear  Band  Development: 

The  presence  of  filter  paper  was  found  lo  hide,  inhibit  and  direct  the  formation  of 
shear  bands.  At  failure,  the  slits  in  the  filter  paper  have  enough  influence  so  that  they 
dictate  the  final  shear  band  pattern. 

The  shear  bands,  whether  they  occur  near  the  peak  in  an  uimotched  specimen,  or  at 
80  percent  of  the  peak  triggered  by  the  presence  of  a  crack,  eventually  take  the  same 
direction  showing  once  more  that  the  behavior  is  dictated  by  the  plasticity  of  the  material. 
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This  direction  was  found  to  be  fixed  with  respect  to  the  direction  of  the  major  principal 
stress.  A  similar  result  was  foimd  by  Arthur  and  Ehinstan  for  sand  materials. 

If  a  specimen  is  subjected  to  further  rotation,  at  or  beyond  the  peak,  an  additional 
shear  band  develops  in  the  vicinity  of  the  crack  tip.  The  first  shear  band,  often  hindered  by 
the  boundary,  stops  and  additional  deformations  become  primarily  localized  in  the  new 
shear  band.  The  second  band  usually  starts  with  a  flat  slope  but  later  tends  to  become 
parallel  to  the  first  one. 

The  inclination  of  the  cracks  was  found  to  act  in  unison  with  the  direction  of 
anisotropy  and  that  of  the  applied  shearing  stress.  The  development  of  the  shear  bands 
seems  to  ignore  the  presence  of  an  upstep  crack  (-ve  a).  Two  bands  will  pass  through  the 
tips  while  others  will  form  in  other  locations;  some  stopping  while  attempting  to  cross  the 
crack.  Stress  concentrations  at  the  tips  of  the  crack  have  little  or  no  influence.  On  the 
other  hand  a  downstep  crack  (>ve  a)  will  become  a  part  of  the  shear  band  and  this 
combination  will  tend  to  weaken  the  specimen.  In  particular,  the  crack  at  22.5  degrees 
blended  so  well  with  the  directions  dictated  by  the  plasticity  of  the  material  that  only  one 
band  was  seen  to  develop  in  this  particular  specimen  to  dissipate  the  energy. 

Isotropic  specimens  give  diffuse  shear  bands.  It  appears  that  the  random 
orientation  of  particles  and  clusters  inhibits  the  formation  of  the  clearly  defined  surfaces 
observed  when  testing  Kq  consolidated  clays.  Also,  the  shear  bands  that  develop  in 
overconsolidated  clays  are  thinner  and  more  pronounced  than  their  counterpart  in 
normally  consolidated  clays;  localization  being  more  intense. 
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Most  of  the  photographs  were  made  at  the  end  of  the  tests  showing  that  parts  of  the 
specimens  between  the  shear  bands  suffer  little  or  no  permanent  deformation.  This  leads 
us  to  question  the  validity  of  calculated  strains  from  measurements  of  deformations  made 
on  the  whole  specimen.  However,  since  shear  bands  occur  at  or  beyond  the  peak  for 
unnotched  specimens,  it  is  only  the  interpretation  of  that  part  of  the  behavior  that  is  in 
question;  for  notched  specimens,  block  motions  start  occurring  at  80  percent  of  the  peak’s 
value. 

While  there  is  a  substantial  amount  of  rebound  that  takes  place  upon  removal  of  the 
stresses,  some  rebound  was  observed  to  take  place  in  the  vicinity  and  upon  formation  of  the 
shear  bands.  The  model  of  Rice  and  Rudnicki  (1980)  takes  into  account  this  possibility. 

Crack  propagation,  shear  banding  and  damage  zones  occur  under  symmetric  cyclic 
loading  for  both  normally  consolidated  and  overconsolidated  clays.  The  bands,  in  a 
herringbone  pattern,  form  a  damage  zone  on  both  sides  and  in  front  of  a  propagating  crack. 
Isotropy  tends  to  narrow  this  zone  and  diffuse  it  such  that  the  bands  are  not  clearly 
delineated.  The  same  happens  with  inclined  cracks  in  anisotropic  specimens;  showing  the 
importance  of  fabric,  direction  of  stress  and  direction  of  cracks. 

Static  and  slow  cyclic  loading  tests  led  to  the  formation  of  shear  bands  and  zones 
where  most  of  the  deformation  at  failure  took  place,  leaving  the  rest  of  the  specimen 
"relatively"  undisturbed.  This  is  not  the  case  for  impulse  loading  where  severe  distortions 
occurred  in  all  the  specimen  and  not  only  in  the  vicinity  of  the  crack.  Viscous  effects  are 
indeed  quite  important  and  must  be  taken  into  consideration  when  modelling  shear  band 
formations  in  saturated  clays. 
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When  two  cracks  are  present,  they  act  independently  or  interact  depending  on  their 
relative  position  and  the  direction  of  the  applied  stresses.  In  our  case  cracks  in  the 
downstep  position  acted  as  a  single  unit.  The  space  between  them  could  suffer  instabilities 
without  affecting  the  overall  behavior  of  the  specimen. 

Since  shear  bands  were  found  to  propagate  in  a  stable  manner,  measurements  can 
be,  and  have  been  made  along  the  crack  and  along  the  shear  bands.  By  extrapolation,  one 
can  find  with  reasonable  accuracy  the  point  at  which  banding  begins  in  the  vicinity  of  a 
crack;  its  evolution  can  thus  be  studied.  Such  measurements  are  indispensable  in  the 
formulation  of  a  bifurcation  and  localization  model. 

7J  In  the  Area  of  Modelling: 

It  was  found  that  models  from  linear  fracture  mechanics  were  not  applicable 
because  of  the  plasticity  of  the  material.  Linear  elastic  analyses,  however,  give  good 
indications  about  the  zones  of  high  tensile  stresses;  but  an  elasto-plastic  analysis  was 
needed  to  explain  the  direction  of  the  shear  bands.  The  Drucker-Prager  model  was 
satisfactory  in  predicting  the  direction  of  the  observed  shear  bands,  proving  once  more  that, 
even  when  tested  undrained,  saturated  clays  behave  as  frictional  materials.  There  is  no 
doubt  that  along  the  shear  bands  the  dilation  angle  is  not  negligible  and  that  water 
migration  takes  place  in  this  zone.  It  is  expected  that,  soon,  local  measurements  of  water 
content  in  the  cell  will  be  possible.  A  combination  of  measurements  and  finite  elements 
analyses  will  be  needed  to  investigate  the  physical  and  mechanical  properties  in  the  vicinity 
of  shear  bands  and  cracks. 
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7.4  In  the  Area  of  Fabric; 


The  measurements  made  with  the  surface  analyzer  and  with  the  Electron 
Microscope  indicate  that  the  deformations  occur  primarily  at  the  level  of  the  clusters.  The 
shear  band’s  thickness  for  EPK  clay  was  estimated  to  be  about  0.5mm;  and  that  is  about 
one  order  of  magnitude  larger  than  the  size  of  a  cluster. 

In  conclusion,  it  is  believed  that  this  study  has  opened  an  important  area  of 
investigation  for  clay  soils.  It  has  shown  that,  with  very  good  but  not  necessarily  very 
expensive  equipment  one  can  go  a  long  way  in  obtaining  the  information  needed  to 
characterize  and  model  materials  behavior. 
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